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optimized.  A simplified  chemical  etching  technique  was  established  for  the 
manufacture  of  such  fluidic  elements. 
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SUMMARY 


A detailed  experimental  and  theoretical  study  of  the 
flow  field  in  the  interaction  region  of  a fluidic  wal 1 -attachment 
amplifier  was  made.  Velocity  profiles  were  measured  using  a 
Laser  Doppler  Velocimeter.  A mathematical  model  of  the  flow  was 
also  established,  using  an  integral  method.  Reasonable  agreement 
between  theory  and  experiments  was  obtained. 

As  a continuation  of  previous  work  (Ref. 11),  the  design 
of  two  bistable  amplifiers,  operating  respectively  with  turbulent 
and  with  laminar  flow,  was  optimized.  A simplified  chemical  etching 
technique  was  established  for  the  manufacture  of  such  fluidic 
el ements . 
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PART  I 


STUDY  OF  THE  JET-TO-WALL  ATTACHMENT 
IN  BISTABLE  AMPLIFIERS 

1.  INTRODUCTION 

Since  the  beginning  of  the  fluidic  technology,  the 
wall  attachment  bistable  amplifier  has  been  one  of  the  units 
most  widely  used  in  fluidic  logic  circuits.  It  has  also  been 
the  object  of  numerous  experimental  and  theoretical  investigations. 
The  recent  references  1 and  2 review  the  state  of  the  art  in 
relating  the  static  and  dynamic  performances  of  the  bistable 
amplifier  to  its  geometry  and  operating  conditions.  Reference 
2 indicates  that  a satisfactory  prediction  of  the  overall 
performances  of  the  bistable  amplifier  can  be  achieved  using  the 
theory  proposed  by  Goto  and  Drzewiecki  (Ref. 3). 

What  is  lacking,  however,  is  a detailed  and  accurate 
description  of  the  flow  field  inside  the  fluidic  element.  Such 
flow  field  is  characterized  by  the  presence  of  a curved  turbulent 
jet  moving  towards  a side  wall  and  attaching  to  it  by  Coanda 
effect  (Fig.l).  Existing  theoretical  descriptions  are  based 
(Ref.l)  on  the  assumption  that  the  jet  develops  as  a free  straight 
jet,  and  therefore  that  its  non-dimensional  velocity  prof i 1 e, spread 
and  centerline  velocity  decay  are  known.  Overall  mass  and  momentum 
balances  are  then  written,  and  the  resulting  system  of  algebraic 
equation  is  solved  for  the  attachment  point,  the  jet  radius  of 
curvature  and  the  separation  bubble  pressure.  However,  in  all  the? 
theoretical  treatments,  some  of  the  following  simplifying  hypotheses 
are  introduced  : the  jet  curvature  and  the  bubble  pressure  are 
constant  in  the  entire  flow  field,  only  the  attachment  wall  is 
present;  the  velocity  and  relative  pressure  in  the  region  outside 
the  convex  boundary  of  the  jet  are  zero;  physically  unacceptable 
flow  patterns  are  thus  often  assumed  or  implied. 

The  need  for  a theory  which  is  not  based  on  those 
simplifying  hypotheses  was  thus  felt,  and  originated  the  work 
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here  reported.  In  addition,  detailed  Laser-Ooppler  velocity 
measurements  in  the  entire  region  between  the  two  attachment 
walls  were  needed  to  provide  a comparison  with  such  theory. 
Pressure  distributions  on  the  attachment  walls,  also  required 
for  the  same  purpose,  had  already  been  measured  in  Ref. 4, 
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2.  EXPERIMENTAL  STUDY 

2 . 1 Research  programme 

The  object  of  the  study  was  to  determine  the  steady- 
state  flow  field  in  the  interaction  region  (Fig.l)  of  a bistable 
amplifier  model.  This  was  done  by  measuring,  along  several 
sections,  perpendicular  to  the  attachment  wall,  the  velocity 
component  parallel  to  the  same  wall.  In  order  not  to  perturbate 
the  flow,  a non-obtrusi ve  instrumentation,  the  Laser  Doppler 
Velocimeter,  had  to  be  used  for  such  measurements.  As  a first 
approach,  the  investigation  was  conducted  on  bistable  configurations 
without  vents  and  splitter,  operated  with  no  control  flow  injection. 
Two  modes  of  operation  were  used  with  control  channels  either 
closed  (Fig. la)  or  open  (Fig. lb)  to  atmosphere. 

2 . 2 Bistable  model 

One  of  the  models  used  in  Ref. 4 was  again  employed 
for  the  present  research.  It  is  schematically  shown  in  Fig. 2. 

In  order  to  get  experimental  data  comparable  with  the  prediction 
of  the  theory  (§  1-3)  which  assumes  two-dimensional  flow  and 
very  long  attachment  walls,  it  was  decided  that  the  model  should 
have  a highaspect  ratio  f/W^  to  decrease  the  importance  of  three- 
dimensional  secondary  flows,  and  an  attachment  wall  length  L^/W^ 
significantly  larger  than  the  expected  jet  attachment  distance 
(L,/W.  about  5 to  10).  A model  having  T/W,  = 8 and  L /W  = 17 
was  therefore  selected.  All  the  other  relevant  dimensions  are 
given  in  Fig. 2. 

The  model  was  constructed  jy  sandwiching  aluminum 
blocks  between  two  side  plates.  The  latter  were  made  out  of 
plexiglass,  except  in  the  zone  where  LDV  measurements  had  to  be 
made  and  where  thin  (2mm)  glass  plates  had  to  be  used.  The  aluminum 
blocks  forming  the  attachment  walls  were  interchangeable,  so 
that  the  wall  angle  a and  set  back  A could  be  varied. 


2 . 3 Test  set  up  and  procedure 


For  conducting  the  tests,  a Laser  Doppler  Velocimeter 
developed  at  the  VKI  was  employed.  This  LDV  uses  a 15  mW  Spectra 
Physics  Helium  Neon  Laser,  a beam  splitting  and  focussing  optical 
unit,  a light  collecting  optical  unit  incorporating  a photo- 
detector, and  an  electronic  data  processing  unit.  This  last  item 
was  completely  designed  and  constructed  at  VKI.  It  differs 
substantially  from  commercially  available  units  and  is  described 
in  references  5 and  6.  Also  the  two  optical  units  were  assembled 
at  VKI  (Ref .5)  . 

For  the  present  research,  the  VKI  Laser  Doppler  Velo- 
cimeter was  operated  in  the  dual-beam,  forward  scattering  mode. 
The  two  beams  were  .40  mm  in. diameter,  and  were  crossing  each 
other  at  an  angle  of  2,73®.  The  Tirobe  volume  was  thus  of 
.40  mm  X .40  mm  x 8.47mm,  the  last  dimension  being  measured  in 
the  direction  of  the  bisectrix  of  the  smaller  angle  formed  by 
the  axes  of  the  two  laser  beams.  Within  the  probe  volume,  the 
interference  fringe  spacing  was  of  6,67  pm.  The  bistable  model 
was  mounted  so  that  the  largest  axis  of  the  probe  volume  was 
perpendicular  to  the  glass  side  plates,  and  was  located  at  mid 
distance  between  them.  The  fringe  planes  were  perpendicular  to 
the  attachment  wall. 

The  air  injected  in  the  supply  nozzle  of  the  model  had 
to  be  seeded  with  small  particles.  These  were  obtained  using  a 
standard  smoke  generator,  producing  a mist  of  condensed  oil 
droplets. 


The  tests  were  conducted  by  displacing  the  LDV  probe 
volume  perpendicularly  to  the  attachment  wall  of  the  bistable, 
and  by  continuously  recording  on  an  X-Y  plotter  velocity  profile. 
Such  measurement  was  repeated  at  various  distances  L measured 
along  the  attachment  wall,  as  shown  in  Fig. 2. 
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Gauge  pressure  in  the  settling  chamber  of  the 
supply  nozzle  was  set  at  100  mm  WG  (981  N/m^).  This  corresponds 
to  a supply  jet  Reynolds  number  of  6660. 

2 . 4 Experimental  results 

The  velocity  profiles,  measured  for  seven  different 
configurations,  obtained  by  varying  the  attachment  wall  angie 
(10®  and  15®),  the  set  back  (0.1;  1 and  2)  and  the  opening  or 
closing  of  the  control  channels  are  given  by  the  solid  line 
curves  in  figures  3 to  9. 

In  such  figures,  u is  the  velocity  component  parallel 
to  the  attachment  wall,  u^  = / 2 p^/p  is  a reference  velocity, 
the  coordinate  y = y"  is  the  normal  distance  from  the  attachment 
wall.  The  distance  L measured  along  the  attachment  wall,  is 
expressed  in  terms  of  the  supply  nozzle  width  . 

Because  the  LDV  used  does  not  detect  the  direction 
of  the  velocity,  but  only  its  magnitude,  the  recorded  profiles 
do  not  show  the  reverse  flow  existing  on  both  sides  of  the  supply 
jet.  In  those  two  regions,  the  true  trend  of  the  velocity  has  been 
obtained  by  manually  drawing  the  dotted  line  curves,  symmetrically 
to  the  recorded  (solid)  ones. 

The  dash  and  dot  curves  appearing  in  the  same  figures 
will  be  described  later  (§  1.3). 

Various  velocity  profiles  have  been  assembled  in  Fig. 10 
to  give  a complete  flow  picture. 

The  data  shows  the  existence  in  the  jet,  close  to  the 
nozzle  exit,  of  a potential  core,  where  the  velocity  is  approximately 
constant.  This  core  is  bounded  by  two  shear  layers  which  are 
themselves  adjacent  to  reverse  flow  regions.  The  potential  core 
extends  downstream  of  the  nozzle  exit  section  for  a distance  of 


3 or  4 nozzle  widths.  Further  downstream  of  this  point  , the  two 
shear  layers  join  together  along  thejet  centerline.  This  latter 
is  defined  as  the  locus  of  points  where  the  velocity  has  a maximum 
amplitude  or,  inside  the  potential  core  as  the  line  dividing  the 
core  into  two  parts  of  equal  width. 

Furthermore,  the  width  of  the  reversed  flow  region 
existing  between  the  jet  and  the  attachment  wall  decreases  to 
zero  when  moving  from  the  nozzle  exit  to  the  attachment  point. 

From  each  velocity  profile,  the  jet  centerline  velocity 
and  the  jet  width  can  be  derived.  If  similarity  exists,  these 
quantities  are  sufficient  to  define  the  velocity  profile.  The 
jet  width  is  usually  defined  as  the  distance  between  the  jet 
centerline  and  the  point  where  the  velocity  reaches  a specified 
value,  for  instance,  that  of  the  reverse  flow  (outside  the  jet), 
plus  0.5  or  0.1  times  the  maximum  velocity  difference  (centerline 
minus  outside).  This  can  of  course  be  done  on  each  side  of  the 
jet;  thus,  four  parameters  b^.  g,  b^.  b^  b^  ^ can  be  defined. 

The  jet  centerline  velocity  and  widths  were  thus 
derived  from  the  measured  data;  they  are  given  as  a function  of 
the  downstream  distance  in  figures  11  to  12  and  13  to  18 
respectively  (open  symbols). 

For  a free  jet,  the  parameters  b ^ and  b ^ vary  linearly 
with  X/W^.  This  is  shown  in  figures  13  to  18  by  the  dotted  straight 
lines.  Comparison  with  the  present  data  shows  that  all  along  the 
outer  side  as  well  as  along  the  initial  part  of  the  inner  side 
of  the  attaching  jet,  the  spread  is  the  same  as  that  of  a free 
jet. 


3.  THEORY 


3 . 1 Statement  of  the  problem 

The  objective  of  the  theory  is  to  describe  mathematically 
the  two  kinds  of  two-dimensional,  incompressible,  viscous  flows 
depicted  in  figures  la  and  lb.  Such  flows  are  of  the  jet  type,  or 
boundary  layer  type  in  the  sense  that  the  longitudinal  component 
of  the  velocity  is  much  larger  than  the  transversal  component 
and  that  the  transversal  velocity  gradient  is  much  larger  than 
the  longitudinal  one. 

In  the  present  case,  the  flow  is  characteri zed  by  the 
existence  of  curved  turbulent  shear  layers,  with  non  negligeable 
longitudinal  and  transversal  pressure  gradients;  of  a separated 
flow  bubble  with  reverse  flow;  and  of  two  irrotational  flow  regions, 
i.e.  the  potential  core  of  the  supply  jet  and  the  region  beyond  the 
convex  jet  boundary. 

3 . 2 Approach  to  the  problem 

The  mathematical  model  here  described  is  based  on  an 
integral  method,  which  assumes  a realistic  shape  for  the  transversal 
profile  of  the  longitudinal  component  of  the  velocity.  This  velocity 
profile,  defined  in  terms  of  unknown  parameters  (widths  and 
maximum  velocities)  is  used  to  integrate  the  mass  and  momentum 
conservation  equations.  The  latter  are  written  by  using  Prandtl 
constant  diffusivity  theory  to  compute  the  Reynolds  stresses, 
and  by  neglecting  the  variation  across  the  jet  of  the  streamlines 
curvature.  Additional  inviscid  flow  relationships  are  written  for 
the  two  irrotational  flow  regions.  The  rate  of  increase  of  the 
two  jet  half-widths  is  assumed  to  be  the  same  as  for  a free  jet. 

Two  geometrical  relations  for  the  widths  of  the  two  reverse  flow 
regions  existing  on  each  side  of  the  jet  are  also  written. 

The  resulting  system  of  equations  is  then  solved  for 
the  velocity  profile  parameters,  as  well  as  for  the  wall  pressure 
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distributions,  the  trajectory  of  the  jet  centerline  and  its 

curvature,  the  attachment  point.  This  is  done  numerically,  by 

dividing  the  flow  field  into  sections  perpendicular  to  the  jet  j 

centerline,  and  by  successively  solving  the  system  of  equations  j 

for  each  section,  assuming  the  curvature  of  the  jet  centerline  j 

and  the  bubble  pressure  to  be  constant  between  each  section 

and  its  neighbouring.  The  calculation  starts  at  the  nozzle  exit 

plane,  and  finishes  in  a section  close  to  the  attachment  point. 

For  the  case  of  closed  control  ports,  the  value  of 
the  bubble  pressure  at  the  nozzle  exit  plane,  not  predicted  by 
the  theory  is  derived  from  experimental  results  (Ref. 4)  and  used 
to  initiate  the  calculation. 

For  the  case  of  open  control  ports,  the  bubble  pressure 
at  the  nozzle  exit  plane  is  given  by  an  additional  expression 
relating  the  pressure  drop  through  the  control  ports  to  the  flow 
entrained  through  the  same  ports.  Continuity  equation  is  also 
used  to  relate  this  flow  to  the  flow  entrained  by  the  jet  and 
to  the  return  flow  existing  out  of  the  jet.  Therefore  the  pressure 
at  nozzle  exit,  on  both  sides  of  the  jet,  as  well  as  the  flows  ! 

I 

entrained  through  both  control  ports  are  also  predicted  by  the  \ 

theory. 

I 

i 

J 

3 . 3 Choice  of  the  velocity  profile 

A curvilinear  orthogonal  coordinate  system  is  used. 

The  X-coordinate  is  measured  along  the  curved  jet  centerline  1 

starting  from  the  nozzle  exit  plane  and  the  Y-coordinate  i 

perpendicularly  to  it. 

The  profile,  in  the  Y-direction,  of  the  X-component 
of  the  velocity  is  assumed  to  be  that  shown  in  figure  19.  This  ‘ < 

is  justified  by  the  experimental  results  shown  for  instance  in 
figure  10. 

. 


The  prof i le  i 5 made  of f i ve  layers  : two  reverse  flow 
regions  with  uniform  velocity,  two  shear  layers,  where  the 
velocity  profile  is  supposed  to  be  similar  to  that  of  a straight 
jet  in  a counterflowing  stream  (Ref. 7),  and  a potential  core 
region,  existing  only  close  to  the  nozzle  exit  plane  where  the 
flow  is  irrotational  and  the  velocity  profile  almost  flat. 

The  velocity  profile  is  therefore  determined  by  9 
parameters,  i.e.  5 character! sti c widths  and  4 characteristic 
velocities. 

The  law  describing  the  velocity  variation  in  the 
potential  core,  i.e.  for  - b^  < Y < b^,  is  obtained  by  integrating 
the  inviscid  Euler  equation  in  the  Y-direction,  i.e.  3p/3y  = pU2/R. 
Then,  assuming  the  radius  of  curvature  R of  the  jet  streamlines 
to  be  constant  along  y,  one  gets  : 


u = U . e 

max  ,e 


(b-y)/R 


-(b  +y)/R 

u = U . . e ^ 
max . 1 


and  therefore 


U . = U . e 

ma  X . 1 max . e 


2b,/R 


Of  course  in  the  established  flow  region,  where  b^  has  decreased 


to  zero,  one  gets  U ^ = U = U 

max.1  max.e  r 


In  the  inner  (i)  or  outer  (e)  shear  layer,  the  velocity 
variation,  suitably  non-dimensional i zed  as 


10 


1 


u - U 


U 


= ^(n).  where  n = — — (k  = i,e)  (3a) 


max . k r .m . k 


win  be  expressed  by  the  relation,  proposed  by  Abramovich  (Ref. 7) 
2 


f(n)  = (1  - 


(3b) 


This  expression  approximates  reasonably  well  (Fig. 20), 
the  exact  solution,  due  to  Goertler,  to  the  boundary  layer  equations 
written  for  a free,  turbulent,  two-dimensional  jet  (Ref. 8). 
Furthermore,  it  does  not  involve,  as  does  the  Goetler  solution 
hyperbolic  functions,  and  therefore  leads  to  much  simpler 
calculations. 


Finally,  it  has  to  be  stressed  that  the  two  velocities 

U and  U . of  the  reversed  flow  outside  the  jet  are  obviously 
rme  rmi 

always  negative  in  the  case  of  closed  controls,  while  in  the  case 
of  open  controls,  can  be  positive,  in  the  region  close  to 

the  inner  control  part.  This  is  due  to  the  fact  that  ambient 
air  is  entrained  into  the  separation  bubble  through  this  inner 
control  part. 


3 . 4 Summary  of  unknowns 

As  previously  explained,  the  flow  equations  can  be 
integrated  once  the  velocity  profile  has  been  defined.  This  leads 
for  the  closed  control  case  to  a system  of  equations  relating 
the  following  unknowns  : 


5 velocity  profile  widths 


U. 


U 


4 velocities 
2 pressures  . 

1 radius  of  curvature  R. 


, u 


rmi  ' max.1  max.e  r.me 
p.  (in  the  separation  bubble),  p 


(outside  the  jet) 


In  the  case  of  open  controls,  two  additional  unknown 
quantities  are  introduced  i.e.,  the  volume  flow  rates 


which  are  entrained  through  the  control  ports. 


The  independent  variable  is  the  abscissa  measured 
along  the  jet  centerline,  or  alternatively,  the  angle  4>  between 
the  nozzle  exit  plane  (x=0)  and  a flow  section  plane,  which  is 
perpendi cul ar  to  the  jet  centerline.  Also  the  distance  along  the 
attachment  wall  L,  will  be  used  as  an  independent  variable. 

3 . 5 Definition  of  the  equations 

3.5.1  Ihe_cl osed_control _case 

The  12  equations  relating  the  12  unknowns  were  derived 
in  Ref. 9.  One  of  them  has  already  been  written  (Eq.2).  The  other 
ones  are  listed  in  appendix  A-1.  Their  physical  meaning  is  given 
here  below  : 

Eq.(4)  : the  width  b^  is  proportional  to  X,  as  for  a free  jet 

Eq.(5)  : the  width  b^  is  proportional  to  X,  as  for  a free  jet, 
if  > 0;  otherwise  this  equation  is  replaced  by  the 
identity  = 0 

Eq.(6)  : the  width  is  given  by  the  available  space  between 
the  outer  jet  edge  and  the  wall  (Fig. 21) 

Eq.(7)  : the  width  (if  > 0)  or  b^.  (if  = 0)  is  given  by  the 

available  space  between  the  inner  jet  edge  and  the  attach- 
ment wall  ( Fi g . 2 1 ) 

Eq.(8)  : the  velocities  U and  U are  related  by  a mass 

' ' max.e  rme  

conservation  condition  equating  the  reverse  flow  existing 
outside  the  jet  to  the  flow  entrained  into  the  jet  along 
its  outer  edge 

Eq.(9)  : the  velocities  U . and  U are  related  by  a mass 
' ' max.i  rme  

conservation  condition  similar  to  the  previous  one,  and 
written  for  the  inner  side  of  the  jet. 
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Eq.(lO) 


Eq.(ll) 


Eq.(12) 


Eq.(13) 


Eq.(14) 


Bernoul 1 i ' s equation  is  used  to  relate  the  pressure 

and  the  velocity  U , in  the  irrotational  reversed 

flow  region  existing  outside  the  convex  jet  boundary 

An  additional  relation  between  L)  and  U is 

— max.e  rme 

obtained  by  combining  the  previous  equation  with 
Bernoul 1 i ' s equation  written  at  the  outer  edge  of  the 
jet  potential  core,  and 

a momentum  balance  across  the  outer  shear  layer  of  the  jet 
The  bubble  pressure  p^  is  obtained  from  a momentum 
balance  across  the  jet,  from  the  dividing  streamline 
to  the  outer  jet  edge.  The  dividing  streamline  is 
defined  by  a mass  conservation  condition  equating,  in  the 
region  bounded  by  the  jet  centerline  and  the  attachment 
wall,  the  jet  forward  flow  to  the  nozzle  supply  flow. 

In  this  case  of  closed  control  ports,  the  dividing 
streamline  coincides  with  the  attachment  streamline, 
bounding  the  separation  bubble 

A condition  of  conservation  of  momentum  along  the 
X-di recti  on  is  written  for  a control  volume  defined  by 
the  nozzle  exit  plane,  the  jet  boundaries,  and  an 
arbitrary  cross  section  of  the  flow 

A momentum  balance  along  the  direction  of  the  attachment 
wall  is  written  for  the  flow  in  the  recirculating 
bubble , bounded  by  two  adjacent  cross  sections  of  the 
flow. 


It  is  important  to  note  that,  for  the  first  flow  section 
X = 0,  the  system  of  12  equations  reduces  to  7 identities,  each 
specifying  the  value  of  one  unknown  and  to  4 equations  relating 
the  remaining  5 unknowns.  Indeed  Eq.(14)  loses  its  meaning  and 
disappears  from  the  system.  Therefore,  the  value  of  one  of  the 
unknowns,  p^ , has  to  be  arbitrarily  specified  to  start  the 
calculation  at  section  X = 0,  and  the  results  obtained  are  then 
used  for  the  calculation  at  the  next  section  and  the  following 


ones . 


I m wm  ■ 
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3.5.2  Ihe_OBen_contrgl _case 

The  main  difference  between  this  case  and  the  previous 

one  is  that  two  additional  unknowns  have  to  be  added,  and 

, i.e.  the  rate  at  which  the  fluid  is  entrained  through  each 

control  port.  As  a result,  the  problem  has  to  be  solved  in  two 

steps  : first  a system  of  14  equations  in  the  14  unknowns,  written 

for  an  initial  control  volume  is  solved,  thus  yielding  in  particular, 

the  values  of  Q and  Q ..  Then  the  second  step  is  involving 
ce  Cl 

the  solution,  at  successive  flow  sections,  of  a system  of  12 
equations  in  12  unknowns,  as  in  the  closed  control  case. 

Furthermore,  in  this  case  the  dividing  streamline 

and  the  attachment  streamline  do  not  coincide.  The  separation 

bubble  is  bounded  by  the  attachment  streamline.  The  pressure  in 

the  bubble  is  assumed  to  be  the  same  as  the  pressure  along  the 

attachment  streamline.  The  flow  scheme  shown  in  Fig. 22  is  adopted. 

As  it  may  be  seen,  the  separation  bubble  does  not  begin  at  the 

control  port  edge,  as  in  reality,  but  at  the  flow  section  where 

the  calculated  value  of  the  velocity  H . is  zero. 

rmi 

3. 5.2.1  System  of  equations  for  the  initial  control 
vol ume 

The  initial  control  volume  is  defined  in  Fig. 22.  It  is 
bounded  by  the  nozzle  exit  plane  (X  = 0)  and  by  that  flow  section 
intersecting  the  attachment  wall  at  the  control  port  edge. 

Nine  of  the  fourteen  flow  equations  are  identical 
respectively  to  equations  ( 2 ) ( 4 ) ( 5 ) ( 6 ) ( 7 ) ( 1 0 ) ( 1 1 ) ( 1 2 ' ) ( 1 3 ) , 
written  for  the  closed  control  case.  The  remaining  five,  derived 
in  Ref. (9),  are  listed  in  Appendix  A-2.  Their  physical  meaning 
is  gi ven  here  below  ; 
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Eq.(15)  : A momentum  balance  along  the  direction  of  the 

attachment  wall  is  written  for  the  flow  bounded  by  the 
dividing  streamline  and  the  attachment  wall. 

This  equation  corresponds  to  Eq.(14)  written  for  the 
closed  controls  case 

Eq.(16)  : A mass  conservation  condition  equates  the  flow  entrained 
into  the  jet  along  its  outer  edge  to  the  sum  of  the 
control  fow  and  the  reverse  flow  existing  outside  the 
jet.  This  equation  corresponds  to  Eq.(8)  written  for 
the  closed  controls  case 

Eq.(17)  : A mass  conservation  condition  similar  to  the  previous 
one,  is  written  for  the  inner  side  of  the  jet.  This 
equation  corresponds  to  Eq.(9)  written  for  the  closed 
controls  case 

Eq.(18)  : The  static  pressure  outside  the  curved  outer  boundary 
of  the  jet,  p , is  given  by  Bernoulli's  equation, 
taking  into  account  the  total  pressure  loss  in  the  outer 
control  channel  . This  pressure  loss  is  a function  of 
the  rate  at  which  the  fluid  is  entrained  through 

the  control  channel.  Such  function  is  obtained  from 
experimental  data 

Eq.(19)  : An  expression  similar  to  the  previous  one  is  written 
for  the  inner  control  channel  to  give  pressure  p^ 

3. 5.2.2  System  of  equations  for  each  successive  flow 
section 

This  problem  is  quite  similar  to  that  described  in 

paragraph  3.5.1.  Indeed,  ten  out  of  the  twelve  equations  have 

already  been  established  (equations  2,4,5,6,7,10,11,13,16,17). 

The  remaining  two  are  ; 


Eq.(16) 


Eq.(17) 


Eq.(18) 


Eq.(19) 


Eq.(21)  ; Same  physical  meaning  as  Eq.(15) 

Eq.(2C)  : The  bubble  pressure  p^.  is  obtained  from  a momentum 

balance  across  the  jet,  from  the  attachment  streamline 
to  the  outer  edge.  The  attachment  streamline  is  defined 
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by  a mass  conservation  condition  equating,  in  the 
region  bounded  by  the  jet  centerline  and  the  attachment 
wall,  the  jet  forward  flow  to  the  nozzle  supply  flow; 
this  equation  corresponds  to  Eq.l2  written  for  the 
closed  control  case 

3 . 6 Presentation  and  discussion  of  results 

A numerical  procedure  was  used  to  solve  the  two 
systems  of  equations  established  in  the  previous  paragraph  for 
the  two  cases  of  closed  or  open  control  parts.  This  was  done  first 
for  the  initial  flow  section  or  the  initial  control  volume,  and 
then  for  each  successive  flow  section. 

Two  computer  programs  were  written,  using  Fortran  IV 
computer  language.  They  are  described  in  Ref. 9.  Their  flow  charts 
and  listings  are  given  in  Appendix  B of  the  present  report.  The 
programs  were  run  on  a CII  Mitra  15  computer,  for  the  7 geometrical 
configurations  of  the  bistable  amplifier  that  had  previously  been 
experi wenta 1 1 y investigated. 

For  the  calculations,  four  constants  had  to  be  specified 
the  jet  spread  coefficient  E = 0.274 
the  jet  momentum  diffusivity  coefficient 

Ki  = 9,77  (in  the  potential  core  region) 

Ki  = 7,67  (in  the  established  flow  region) 
the  coefficient  expressing  the  effect  of  jet  curvature  on  shear 
stress  K2  = 10 

The  values  of  E and  Kj  are  derived  from  well-known  data  for  free 
turbulent  jets.  The  value  of  K2  was  obtained  by  averaging  the 
data  of  several  investigations  made  on  curved  turbulent  jets. 

As  a result  of  the  calculations,  the  nine  velocity 
profile  parameters  and  the  relative  pressures  on  the  two  bistable 
walls  were  obtained.  This  data  was  used  : 
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- to  draw  the  dimensional  velocity  profiles  (Figs  3 to  9).  The 
calculated  velocity  component  parallel  to  the  jet  centerline  is 
given  versus  the  distance  y'  from  the  attachment  wall  measured 
in  the  direction  perpendicular  to  the  jet  centerline. 

- to  plot  the  decay  with  increasing  distance  from  the  nozzle 
of  the  jet  center  line  velocity  (Figs  11  and  12). 

- to  plot  the  increase,  with  increasing  distance  from  the  nozzle, 
of  two  typical  shear  layer  thicknesses  (b  ^ and  b ^),  for  the 

two  sides  of  the  jet  (Figs  13  to  18). 

- to  plot  wall  pressure  distributions  (Figs  23  to  27). 

On  Figs  3 to  18,  were  already  plotted  the  results  of 
the  LDV  measurements,  taken  at  sections  perpendicular  to  the 
attachment  wall,  while  the  calculations  were  made  at  sections 
perpendicular  to  the  jet  centerline.  Thus  experimental  and 
theoretical  data  are  given  in  two  different  reference  systems. 
Comparison  is  nevertheless  meaningful,  because  the  difference 
between  the  two  reference  systems  is  a rotation  of  a small  angle 

(a  - 4>)  . 


On  Figs  23  to  27  are  also  plotted  results  of  measurements 
of  the  wall  pressure  distribution  previously  carried  out  (Ref. 4). 

As  may  be  seen  from  the  figures,  the  overall  trend  of  the  velocity 
field  is  correctly  predicted  by  the  theory.  The  computed  centerline 
velocities  are  slightly  greater  than  the  measured  ones.  Also  the 
calculated  length  of  the  potential  core  is  about  one  nozzle  width 
greater  than  the  measured  one. 

Figs  13  and  14,  show  very  good  agreement  between  the 
experimental  and  the  theoretical  spread  of  the  jet  on  the  outer 
side  of  the  jet.  This  result  could  of  course  be  anticipated, 
because  Eq.(4)  specifically  assumes  that  the  curved  jet  spread 
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Figs  15  and  16  show  that  good  agreement  is  obtained 
for  the  inner  region  of  the  jet  only  far  before  the  attachment. 

This  could  also  be  expected.  Indeed  the  width  b^  ^ , first 

increases  with  X at  the  same  rate  as  a free  jet,  but  then  decreases 
when  approaching  the  attachment  point.  Again  this  is  implied 
by  the  physical  meaning  of  Eq.(5b)  and  (7). 


Figs  17  and  18  do  not  show  for  the  inner  side  of  the 
jet  good  agreement  between  the  predicted  and  measured  values  of 
the  hal f-vel oci ty  shear  layer  thickness  b^  ^ , when  approaching  ' 

the  attachment  point.  Even  the  trends  are  opposite.  This  disagreement  t 

when  compared  with  the  good  agreement  observed  between  the  full 
widths  indicate  that,  when  approaching  the  attachment  point, 

the  true  shape  of  the  velocity  profiles  differs  from  that  assumed  r 

(Eq.3)  in  the  theory. 

Figs  23  and  24,  related  to  the  closed  control  case, 
show  reasonable  agreement  between  theory  and  experiment  except 
for  the  value  of  the  pressure  peak  associated  with  attachment. 

The  abscissa  of  this  peak,  i.e.  the  attachment  point  location  is  ’/ 

correctly  predicted.  Figs  25  to  27  show  that  in  the  open  control  ! 

case  a better  agreement  is  obtained. 

In  the  case  of  open  controls,  the  theory  also  reasonably 
predicts  the  value  of  the  flow  entrained  through  the  control 
channels,  as  shown  by  the  following  table  : 


10°  1 0.152  0.131 
10°  2 0.174  0.161 
15°  1 0.162  0.144 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 

Detailed  data  on  the  flow  field  inside  the  interaction 
region  of  a bistable  amplifier  was  obtained  by  LDV  measurements. 

This  flow  field  was  reasonably  modelled  by  an  integral 

method . 

Improvement  to  the  theory  and  continuation  of  this 
work  should  be  carried  out,  so  that 

- Eq.(3)  describing  the  velocity  profile,  should  be  written 

f(n)  = (1  - n")' 

so  that  the  value  n = 1.5  adapted  for  the  present  theory  could 
be  relaxed.  Thus  the  shape  of  the  profile  could  change  when 
approaching  reattachment. 

- The  effect  of  forced  control  flow  injection  could  be  taken 
in  account  by  adding  an  extra  pressure  term  in  Eq.(19).  Using 
one  of  the  three  usual  switching  criteria,  the  switching  point 
could  thus  be  predicted. 
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DESIGN  AND  CONSTRUCTION  OF  A RISTABLE  AMPLIFIER 


1.  INTRODUCTION 


VJall  attachment  bistable  amplifiers  operatinn  in  the 
turbulent  flov/  regime  are  commercially  available  since  more  than 
ten  years.  They  are  characterised  by  a pressure  recovery  of  about 
35%,  a pressure  gain  which  is  generally  below  5,  a fan  out 
generally  below  4, 


An  experimental  investigation  carried  out  at  VKI 
(Ref. 11),  on  the  mechanisms  of  bistable  switching  showed  that 
these  performances  could  be  improved.  The  research  was  conducted, 
in  a first  phase,  on  bistable  configurations  without  splitter  and 
without  output  ducts.  Its  purpose  was  to  optiru7e  the  geometry 
of  the  two  straight  attachr.ent  walls.  This  was  done  in  the 
turbulent  and  in  the  laminar  flow  renimcs. 


In  the  turbulent  regime,  the  optimal  geometry  was 
obtained  for  the  usual  value  of  the  aspect  ratio  t/w^  = 2,  with 
a wall  angle  a of  10°,  a wall  length  equal  to  f to  f times  the 
nozzle  width,  and  a wall  setback  a of  1/lP  of  the  nozzle  width. 

For  this  optimum  conf  i nura  ti  on  , the  ratio  of  the  supply  nomentun: 
flux  to  the  minimum  control  momentum  flux  required  for  switcliing 
the  device  was  of  4C.  Furthermore,  the  average  jet  dynamic  pressure 
at  the  end  of  the  attachment  wall  was  about  45%  of  the  supply 
total  pressure. 

From  this  data,  it  is  possible  to  estimate  approximately 
the  performance  of  the  best  bistable  which  would  incorporate  the 
optimized  attachment  walls  with  the  addition  of  the  splitter  and 
output  ducts.  Thus  maximum  achievable  pressure  recovery  would  also  be 
of  45%.  For  this  estimation  the  output  ducts  are  supposed  to  collect 
all  the  jet  flow,  and  the  losses  in  the  vent  renion  are  neglected. 
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Furthermore  the  pressure  gain,  approximately  equal  to  the 
momentum  ratio  times  the  estimated  pressure  recovery,  v.'ould  be 
of  about  20.  These  performances  compare  favourably  v/ith  those 
relative  to  a commercial  device. 

The  tests  in  the  laminar  flow  regime  (Ref. 11)  demonstrated 
the  possibility,  theoretically  predicted  by  Goto  and  Trzewiecki 
(Ref. 3)  of  designing  a fluidic  bistab'e  amplifier  based  on  the 
wall  attachment  of  a laminar  jet  The  optimal  geometry  of  the 
attachment  walls  was  found  to  be  the  same  as  in  the  turbulent  flow 
case.  For  such  configuration  the  maxi n; urn  ratio  of  supply  to 
control  momentum  fluxes  at  switching  was  lying  between  300  and  500. 
It  was  therefore  concluded  that  a laminar  bistable  an.plifier  could 
have  a much  greater  gain  than  a turbulent  one. 

On  tlie  basis  of  such  promisinn  results,  it  was  decided 
to  continue  the  research  described  in  Pef.ll  by  designing  two 
bistable  amplifiers  with  splitter  and  with  output  ducts,  operating 
respectively  in  the  laiinnar  and  the  turbulent  regimes. 

The  objective  was  first  to  optimize  experimentally  the 
geometry  of  the  bistable  amplifiers  using  a large  scale  model  with 
interchangeable  and  moveable  walls,  splitter  and  output  ducts. 

In  a second  phase,  real  scale  amplifiers  would  then  be  constructed 
using  a cher.ical  etching  technique,  and  would  be  tested  for 
evaluating  their  performances. 
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2.  DESIGN  STUDY  OF  A TURBULENT  BISTABLE  AMPLIFIER 

2 . 1 Introduction 

The  purpose  of  this  study  was  to  determine  a good 
design  for  a turbulent  bistable  amplifier.  The  assumption  was 
made  that  the  optimization  of  the  attachment  walls  was  independent 
of  that  of  the  splitter  and  output  ducts. 

The  optimization  of  the  attachment  walls  had  already 
been  carried  out  (Ref. 11)  in  the  case  of  straight  walls.  It  was 
therefore  decided  to  complete  this  research  by  investigating 
the  effects  of  wall  curvature.  Thus  the  best  geometry  of  attachment 
walls,  straight  or  curved,  would  be  determined. 

A splitter  and  two  additional  walls  boundinn  the  output 
duets,  would  then  be  added  to  this  best  configuration.  The 
divergence  angle  of  the  splitter  and  of  the  output  ducts,  as 
well  as  the  width  and  offset  of  vents  would  be  varied  to  get 
the  best  design. 

2 . 2 The  experimental  model  and  test  set  up 

Out  of  the  models  used  for  the  previous  experimental 
research  (Ref. 11),  the  one  having  an  aspect  ratio  t/w^  = 2 was 
selected  for  the  present  study.  This  choice  was  made  for  two 
reasons  : first  it  is  a good  bistable  configuration  accordino  to 
Ref. 11;  second  this  aspect  ratio  is  sufficiently  small  to 
guarantee  a reasonable  simplicity  in  the  manufacture  of  a 
miniaturized  amplifier. 

Fig. 28  shows  schematically  the  model,  and  defines  its 
significant  geometrical  parameters.  The  supply  settling  chamber 
and  nozzle,  not  shown  in  Fig. 28,  have  the  same  planform  as  the 
model  shown  in  Fig. 2.  The  control  channels  have  the  same  width 
as  the  supply  channel.  The  splitter  and  the  output  walls  have 
a rounded  leading  edge.  The  attachment  walls,  the  output  walls 
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and  the  splitter  are  moveable,  so  that  a,y  and  can  be  varied. 
The  use  of  interchangeable  pieces  allows  to  vary  also  a , 6 , > 
as  well  as  the  attachment  wall  shape.  Five  wall  shapes,  also  shown 
in  Fig.  28,  were  investigated. 

Wall  type  A has  the  optimal  shape  according  to  Ref. 11. 
The  others  are  derived  from  it  by  : enlarging  and  shaping  the 
vent  duct  (type  B);  rounding  the  leading  edge  (type  F);  replacing 
the  straight  wall  by  a constant  curvature  with  the  leading  and 
the  trailing  edges  at  the  same  location  as  on  the  straight  wall, 
and  having  a leading  edge  angle  of  5°  (type  C);  rotating  wall 
type  C around  its  leading  edge  by  5®  towards  the  centerline  of 
the  bi stabl e ( type  P)  . 

The  main  nozzle  and  one  of  the  control  ports  were 
connected  to  two  independently  regulated  air  supplies.  The  main 
(supply)  and  the  control  flows  were  measured  by  two  rotameters. 

The  output  volume  flow  was  measured  by  using  the 
device  shown  in  Fig. 28,  made  of  a 200mm  long,  20mm  internal 
diameter,  cylindrical  duct  incorporating  a total  head  tube  and 
a static  pressure  hole.  The  resulting  pressure  differential  was 
calibrated  by  comparison  with  a rotameter.  Static  pressures  were 
measured  under  steady  state  conditions  in  the  settling  chamber  of 
the  supply  nozzle,  as  well  as  in  the  control  channel  and  in  one 
of  the  settling  chambers  connected  to  the  output  ducts. 

2 . 3 Tests  with  attachment  walls  only 

For  the  five  types  of  walls  shown  in  Fig. 28,  the  control 
flow  required  to  switch  the  supply  jet  was  measured.  Supply 
pressure  was  varied  from  400  to  4000  N/m^ . This  corresponds  to  a 
supply  Reynolds  number  varying  from  4^200  to  13,340. 

Walls  type  A,  B and  E were  positioned  at  a setback 
(Fig. 28)  equal  to  1/10  of  the  supply  nozzle  width.  For  wall  E, 
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this  led  to  a secondary  throat  width  w'  (Fig. 28)  larger  than 
for  wall  type  A or  B.  For  this  reason  also  a test  with  w'  = w 
was  made,  for  wall  type  E. 

Walls  type  C and  D were  mounted  with  their  leading 
edge  at  the  same  position  as  the  leading  edge  of  wall  type  A. 

Results  are  given  in  Fig. 29  in  terms  of  the  ratio  of 
the  required  control  flow  to  a reference  flow  = tw  / 2p^/p . 

The  latter  is  the  ideal  supply  flow,  i.e.  calculated  assuming 
a nozzle  discharge  coefficient  of  unity. 

The  results  indicate  that  : 

- the  change  from  wall  A,  having  the  optimal  shape  according  to 
Ref.  11,  to  wall  E,  having  a shaped  vent  duct  has,  as  could 

be  expected,  no  influence  on  the  switching.  Furthermore,  when 
increasing  the  supply  Reynolds  number,  the  switching  flow  ratio 
first  varies,  then  rapidly  reaches  a constant  value.  The  latter 
is  of  0.12  for  Re  > 7000  approximately.  Taking  into  account  the 
actual  value  of  the  discharge  coefficient  of  the  supply  nozzle 
(measured  as  = 0.86)  and  an  estimated  value  of  the  output  to 
supply  flow  recovery  ratio  (around  1.5),  the  previous  figure 
leads  to  an  estimated  flow  gain  of  the  bistable  amplifier  of 
about  11. 

- with  curved  walls  like  type  C or  type  D,  the  same  trend  is 
observed.  However  a constant  switching  flow  ratio  is  obtained 
only  for  Re  > 10,000  approximately.  Furthermore,  adopting  such 
kind  of  walls  instead  of  the  previous  ones  does  not  seen; 
advantageous,  because  it  causes,  in  the  region  of  Reynolds 
number  independence,  an  increase  in  the  control  flow  requirement. 

- straight  walls  with  rounded  leading  edge  (type  E),  when  positioned 
with  the  same  setback  as  walls  A,  yield  a switching  flow  ratio  of 
about  0.2  for  Re  > 8000.  This  indicates  a larger  control  flow 
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requirement  than  in  case  A.  When  the  secondary  throat  width,  and 
therefore  the  setback,  are  decreased,  the  required  central  flow 
decreases,  but  remains  dependent  on  Reynolds  number.  This  of  course 
is  an  undesirable  feature. 

One  conclusion  from  these  tests  can  be  the  adoption 
of  straight  walls  (type  A or  B) , as  a compromise  between  high 
flow  gain  and  Reynolds  number  independence. 

An  additional  remark  can  be  made  : in  Ref. 2,  a modified 
Reynol ds  number  Re ' . 

R4  = Rg/(t^/W^  + 1)(1  + w^/t)2 

is  proposed  for  assessing  the  condition  of  Reynolds  number 
i ndependence  , i .e  . 

> 1000 

In  the  present  case  the  supply  nozzle  and  contraction  are 
quite  long.  They  were  designed  in  such  a way  to  yield  a very 
regular  velocity  profile.  But  when  the  previous  criterion  for  the 
minimum  Reynolds  is  applied,  one  gets  = 12  and  t/w^  = 2) 

R„  > 29.3  X 1000 
e 

This  figure  is  in  apparent  conflict  with  some  of  the 
results  shown  in  Fig. 29,  specially  for  the  use  of  straight  walls, 
where  the  condition  Re  > 7000  was  found.  On  the  other  hand,  it 
justifies  the  result  relative  to  wall  E with  w'  = w^ . Furthermore, 
it  suggests  that  nozzles  much  shorter  than  that  shown  in  Fig. 2 
should  be  used  for  the  final  design  of  a bistable  amplifier. 
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2 . 4 Effects  0'^  splitter  and  output  ducts 


Using  the  straight  walls  (type  A)  discussed  in  the 
previous  section,  46  different  bistable  amplifiers  were 
assembled  by  varying  the  vent  width  w , the  output  walls  set- 
back y and  angle  o and  the  splitter  angle  6.  The  leading  edge 
of  the  splitter  was  located  at  the  same  downstream  distance  from 
the  nozzle  as  the  leading  edges  of  the  output  walls. 

Each  configuration  was  tested  to  determine  its  output 
characteristic  curve.  The  latter  relates  the  output  pressure  and 
the  output  volume  flow.  All  points  on  such  curve  are  obtained 
by  progressively  closing  the  active  outlet.  Thus  output  flow 
decreases  and  output  pressure  increases.  Under  certain  circum- 
stances the  maximum  output  pressure  with  zero  output  flow  cannot 
be  obtained,  because  the  jet  switches  to  the  unrestricted  outlet. 

This  undesirable  phenomenon,  known  as  1 oa d- s wi tch i n g , was  found 
to  occur  for  22  bistable  configurations.  For  the  remaining  24 
configu rations,  the  pressure  (Rp)  and  flow  (R^)  recoveries,  de- 
fined as  the  maximum  ratio  of  an  output  quantity  to  the  corres- 
ponding supply  quantity,  were  derived  from  the  measured  curves. 

i 

The  tests  were  conducted  with  control  ports  either 
both  closed  or  both  ooen  to  atmosphere. 

The  results  are  summarized  in  the  five  following 
tables.  They  indicate  which  configurations  are  switching  under 
load  (LS)  and  for  the  others,  give  the  pressure  recovery  (fioures 
on  upper  line)  and  the  flow  recovery  (fioure  on  lower  line).  When 
two  figures  are  given  for  the  pressure  recovery,  the  first  one 
replies  to  the  closed  control  case,  the  second  one  to  the  open 
control  case. 
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Inspection  of  these  tables  reveals  that  ; 

- switching  under  load  occurs  when  the  vent  width  is  too  small, 
as  could  be  expected,  or  when  the  output  wall  setback  y is 
too  large.  In  other  terms,  the  vent  width  required  to  avoid 
load  switching  increases  with  increasing  output  wall  setback. 

The  dotted  line  on  the  tables  separates  the  cases  in  which 
load  switching  does  occur  from  the  others. 

- increasing  the  vent  width  causes  as  could  be  expected  a 
decrease  in  pressure  recovery,  due  to  turbulent  mixing  losses, 
and  an  increase  in  flow  recovery,  due  to  the  increased 
entrainment  of  air  through  the  vents. 

- increasing  the  wall  offset  seems  to  cause  an  increase  both 
in  flow  and  pressure  recoveries. 

- it  appears  that  several  designs  give  good  performances.  They 
are  indicated  by  an  asterisk  in  the  tables,  and  are  character- 
ized by  a pressure  recovery  above  .45  and  by  a flow  recovery 
above  1.50. 

Among  these,  the  design  characterized  by  o = 10°, 

L/w^  = 8,  a/w  = 0.1,  6 = 5.5°,  y = 10°,  w /w^  = 2 and  y/w  = 1.0 
was  selected  for  the  construction  of  a fluidic  bistable  amplifier 
with  smaller  dimensions  than  the  model  used. 

For  this  configuration.  Fig. 30  gives  the  measured  output 
characteristic,  in  the  two  cases  of  controls  open  or  closed,  and 
in  the  two  cases  of  straight  or  shaped  vents. 

2 . 5 Conclusions  and  recommendations 

As  a result  of  the  parametric  study  undertaken,  some 
geometrical  parameters  of  a bistable  amplifier  operating  in  the 
turbulent  regime  were  optimized. 

It  was  found  that  straight  attachment  walls  were 
performing  better  than  curved  ones.  Straight  attachment  walls 
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with  rounded  leading  edge  might  yield  even  better  performances, 
though  dependent  on  Reynolds  number,  thus  on  supply  pressure. 

Ways  of  reducing  this  dependence,  eventually  by  modifying  the 
shape  of  the  supply  nozzle  should  be  investigated. 

It  was  found  that  optimum  width  of  vents  and  optimum 
setback  of  output  ducts  were  respectively  of  about  2 and  0.5  to 
1 times  the  supply  nozzle  width.  Optimum  divergence  of  output 
ducts  was  about  2.5°  to  5°. 

Furthermore  shaping  the  vent  channel  (wall  type  P)  has 
no  effect  on  the  input  characteristic  and  improves  the  output 
characteri Stic. 

As  a conclusion  of  this  study,  the  silhouette  of  a 
turbulent  flow  bistable  amplifier,  to  he  constructed  by  a chemical 
etching  technique  (Section  II-4)  was  designed  and  is  represented 
in  Fig. 31.  It  should  have  a maximum  pressure  and  flow  recoveries 
respectively  of  0.47  and  1.50,  together  with  a maximum  flow  gain 
of  11. 
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3.  DESIGN  STUDY  OF  A LAMINAR  BISTABLE  AMPLIFIER 
3 . 1 Introduction 


It  is  well-known  that  laminar  flow  fluidic  devices  | 

have  several  advantages  over  turbulent  ones  : absence  of  random  | 

fluctuations  in  the  flow,  and  therefore  high  signal  to  noise 
ratio  and  accurately  defined  switching  conditions,  low  power 
consumption,  possibility  of  miniaturization. 

The  only  existing  type  of  bistable  amplifier  operating 
with  laminar  flow  is  that  designed  by  Drzewiecki  (Ref. 12).  Its 
operation  is  not  based  on  the  phenomenon  of  jet-to-wall  attachment, 
but  rather  on  the  instability  of  the  power  jet  of  a laminar  pro- 
portional amplifier  (Ref. 13,  14).  The  device  has  a pressure  gain  l 

of  15,  a fan  out  of  8,  a power  consumption  of  3.5  mU',  the  supply  ’ 

nozzle  having  a width  of  1 mm  and  a thickness  of  0.5mm.  | 

The  purpose  of  the  present  work  was  to  design  another  ^ 

type  of  laminar  flow  bistable  amplifier,  based  on  jet  to  wall  i 

attachment,  a phenomenon  which  occurs  not  only  in  the  turbulent  | 

flow  regime,  but  also  in  the  laminar  one  (Ref. 11). 

^ ,1 

In  Ref. 11,  laminar  jet  to  wall  attachment  was  inves-  -i 

tigated  only  in  the  case  of  bistable  configurations  without 

f 

splitter  and  without  output  ducts.  The  aspect  ratio  t/w^  was  of 
8 and  of  2,  and  the  control  ports  were  closed.  To  make  the 

» 

design  practical,  the  case  of  open  control  ports,  with  a splitter  j 

and  two  additional  walls,  bounding  the  output  ducts  had  to  be 

considered.  An  aspect  ratio  of  2 was  selected  for  this  further  ' ' 
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3 . 2 The  experimental  model  and  test  set  up 

The  model  of  aspect  ratio  2 used  in  Ref. 11  was  employed 
for  the  present  study,  with  the  addition  of  three  extra  pieces 
(splitter  and  two  output  walls).  The  geometry  of  the  model  was 
varied  during  the  research;  it  is  defined  by  the  parameters 
previously  used  (Fig. 28). 

As  done  in  Ref. 11,  the  tests  were  conducted  with  water. 
Indeed  with  air  at  low  Reynolds  numbers,  the  velocities  are  too 
small  for  accurate  measurements.  With  water,  the  velocities  are 
15  times  lower  than  with  air,  for  the  same  Reynolds  number,  but 
flow  visualizations  are  much  easier.  Furthermore  the  use  of 
water  leads  to  a more  accurate  measurement  of  switching  times 
because  they  are  15  times  greater  than  in  air. 

The  model  was  placed  in  the  VKI  water  table,  1 cm 
below  the  free  surface.  The  output  settling  chamber  shown  in 
Fig. 28  was  removed. 

The  main  nozzle  was  connected  via  a rotameter  to  a 
regulated  water  supply.  A dye  was  injected  into  the  settling 
chamber  of  the  supply  nozzle  to  visualize  the  flow.  Steady  state 
observations  of  the  flow  were  thus  made  with  the  control  and 
output  ports  independently  open  or  closed. 

Dynamic  tests  were  also  conducted.  For  this  purpose, 
a second  water  supply  was  connected  via  a rotameter,  to  a two- 
way  electrovalve,  normally  discharging  the  flow  through  one  of 
its  outputs  in  the  water  table.  The  second  output  of  the  electro- 
valve, normally  closed,  was  connected  via  a short  line  of  small 
impedance,  to  one  of  the  control  ports  of  the  bistable  model, 
the  other  control  port  being  closed.  Actuating  the  electro -valve 
resulted  in  an  abrupt  injection  of  water  into  one  of  the  control 
ports,  and  this  caused  the  bistable  to  switch. 
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This  mode  of  operation  is  termed  here  "with  controls 
initially  closed".  A second  mode  was  achieved  by  injecting  into 
the  control  port  connected  to  the  el ectro- val ve  an  additional 
flow,  equal  to  the  flow  that  would  normally  be  sucked  into  the 
control  port  if  it  were  open.  A separate  water  supply  was  used 
for  this  purpose.  Actuating  the  el ectro-val ve  resulted  in  an 
abrupt  increase,  above  an  initial  non-zero  value,  of  the  flow 
rate  injected  into  one  of  the  control  ports.  The  other  control 
port  was  kept  open  during  all  this  process,  and  this  mode  of 
operation  is  termed  here  " with  controls  initially  open". 

The  flow  in  the  fluidic  device  was  visualized  by 
injecting  a dye  in  the  settling  chamber  of  the  supply  nozzle 
and  by  using  water  of  a different  colour  to  feed  the  control 
por t . 

The  switching  of  the  coloured  water  jet  was  detected 
by  a 1 amp-and-photodiode  system,  located  across  the  initially 
inactive  output  channel.  The  switching  time  was  measured,  using 
an  electronic  counter,  as  the  time  elapsing  from  the  activation 
of  the  el ectro- val ve  to  the  detection  of  the  passage  of  the 
coloured  water  in  front  of  the  photo-diode. 

3 . 3 Steady  state  test  results 

Steady  state  laminar  flow  tests  were  conducted  first, 
with  the  purpose  of  determining  a configuration  exhibiting  a 
bistable  behaviour  either  with  controls  open  or  closed,  and  this 
independently  of  the  loading  or  unloading  of  the  outputs. 

Fig. 32  shows  the  geometry  of  the  configuration  investi 
gated  first.  With  a wall  angle  of  15°,  the  setback  a was  varied 
from  0.1  w^  to  1 w^,  and  the  wall  length  from  6 to  14  times 
w^ . The  triangular  splitter  was  located  downstream  of  the  vents. 

While  such  configurations  were  bistable  when  the 
controls  were  closed,  they  were  not  operating  satisfactorily 
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when  the  controls  were  open.  Indeed,  two  cases  were  observed  : 

- the  laminar  jet  was  attaching  to  one  wall,  but  further 
downstream,  was  separating  from  it  and  impinging  on  the 
splitter  leading  edge.  Thus  part  of  the  jet  was  diverted 
towards  the  wrong  output.  This  situation,  depicted  in  Fig. 33, 
was  occurring  for  the  longer  attachment  walls. 

- the  laminar  jet  did  not  attach  to  the  walls,  but  proceeded 
straight  on  and  equally  divided  upon  inningement  of  the  splitter. 
This  symmetrical  situation,  arising  with  the  shorter  walls, 

is  depicted  in  Fi g . 34 . 

In  order  to  facilitate  the  laminar  jet  attachment,  a 
lower  wall  angle  (a  = 10°),  as  well  as  another  shape  for  the 
splitter,  as  shown  in  Fig. 35,  were  used,  but  without  success. 

It  proved  necessary  to  decrease  the  wall  angle  to  7° 
to  get  a configuration  exhibitinci  bistable  behaviour  in  both 
cases  of  closed  or  open  control  ports.  This  configuration, 
however,  was  not  yet  satisfactory,  because  it  was  load  sensitive  : 
when  blocking  the  active  output,  switching  occurred. 

To  eliminate  this  phenomenon  of  1 oad - swi tch i no  , two 
modi fi cations  were  devised  : the  increase  of  the  vent  dimension 
w^ , and  the  decrease  of  the  splitter  distance  P^p- 

The  first  action  proved  to  be  effective,  as  shown  in 
Fig. 36.  A stable  position  with  the  active  output  blocked  is 
obtained,  no  flow  being  spilled  through  the  other  output.  However, 
the  required  vent  dimension  to  achieve  load  insensitivity  is  too 
great  (w^  » 32  w^),  and  probably  implies  a low  pressure  recovery, 
or  a large  sensitivity  to  external  perturbations. 

The  second  action  proved  to  be  uneffective,  as  shown 
in  Fig. 37.  Indeed  blocking  the  active  output  causes  part  of  the 
jet  to  be  deviated  towards  the  other  output,  by  the  effect  of 
the  sharp  splitter.  On  the  other  hand  operation  with  outputs 
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unblocked  was  satisfactory. 

Because  of  this  observed  bad  effect  of  the  sharp 
splitter,  a blunt  one  was  tried.  This  gave  good  results,  as  shown 
in  Fig. 38.  Indeed  the  blunt  leading  edge  of  the  splitter  causes 
a vortex  to  form  on  the  separated  side  of  the  laminar  jet.  Such 
vortex  helps  in  keeping  the  jet  attached  to  the  wall,  when  the 
active  output  is  blocked.  This  phenomenon,  occurring  here  in 
laminar  flow  is  similar  to  the  well-known  one  (Ref. 15)  occurring 
in  conventional  turbulent  flow  bi stables  having  blunt  or  cusped 
spl i tters . 

The  geometry  of  the  final  laminar  flow  device,  bistable 
with  control  ports  closed  or  open,  is  shown  in  Fig. 39. 

Because  of  the  expected  difficulties  in  accurately 
manufacturing  the  splitter  in  a miniaturized  device,  a further 
study  was  made  on  the  effect  of  splitter  shape  and  position.  The 
following  conclusions  were  drawn  : 


- the  splitter  has  to  be  blunt,  without  rounded  edges, 

- the  splitter  width  may  be  increased,  except  in  the  vicinity  of 
its  leading  edge,  where  a nominal  length  of  1 w^,  with  a 

width  of  0.8  w^  has  to  be  provided. 

- the  splitter  distance  may  be  varied  from  12.6  to  16.4  times 


3.4  Dynamic  switching  test  results 


For  the  two  bistable  configurations  similar  to  that 
shown  in  Fig. 39  and  having  D^p/w^  = 12.6  and  16.4,  the  switching 
times  were  measured  aS  explained  in  section  3.2.  They  are  given  in 
Figs  40  to  47,  versus  the  amplitude  of  the  step  in  the  control 
momentum  flux.  Such  curves  are  called  "dynamic  switching  charac- 
teristics". The  four  cases  of  operation  with  controls  and  outputs 
open  or  blocked  were  considered.  In  the  case  o*  operation  with 
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controls  "initially  open"  (Section  3.2),  also  the  momentum 
flux  J . of  the  flow  entrained  into  the  control  port  prior  to 
the  application  of  the  control  step,  is  given  in  the  figures. 

The  tests  were  carried  out  forfbur  values  varying 
from  350  to  1500  of  the  Reynolds  number  based  on  the  average 
supply  velocity  and  the  nozzle  width. 

The  minimum  ratio  of  control  to  supply  momentum  flux 
required  for  switching  can  be  derived  from  figures  40  to  47. 

The  reciprocal  of  this  quantity  is  proportional  to  the  maximum, 
or  steady  state,  gain.  Also  the  switching  tine  for  an  arbitrary 
ratio  of  supply  to  control  momentum  flux,  for  instance  1/20,  can 
be  derived  from  the  figures.  The  effect  on  these  two  parameters, 
of  a variation  in  splitter  position  D,  in  supply  Reynolds  number, 
and  in  the  operating  condition  (controls  and  outputs  open  or 
closed)  is  shown  in  the  following  table.  In  the  latter  the  first 
figure  is  1000  ( second  is  ( j ^ = 200  c ) ^ ^ ^ ’ 

Inspection  of  this  table,  and  of  Figs  40  to  47  shows 

that  : 

- the  dynamic  switching  characteristics  are  Reynolds  dependent. 

An  increase  in  Reynolds  nun.ber  causes  a decrease  in  switching 
time  and  in  increase  in  the  minimum  control  to  supply  momentum 
ratio,  and  therefore  a decrease  in  maximum  gain. 

The  dependence  of  the  maximum  gain  on  Reynolds  number  is 
equally  strong  in  all  the  32  cases  investigated  : by  increasing 
the  Reynolds  number  by  less  than  one  order  of  magnitude,  the  gain 
decreases  of  more  than  one  order  of  magnitude.  Only  the  case 
where  D/w^  = 16.4,  Re  = 650  to  900,  with  controls  open,  is  an 
exception  . 

The  dependence  of  switching  time  on  the  Reynolds  number 
is  stronger  when  the  splitter  is  located  closer  to  the  nozzle. 
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Operating  the  device  with  the  controls  initially  open  or  closed 
(i.e.  respectively  with  or  without  vents  in  the  control  channels) 
does  not  modify  this  Reynolds  number  dependence.  On  the  other 
hand,  blocking  the  outputs  causes  a decrease  in  this  dependence  in  ■ 
the  case  of  the  shorter  splitter  distance  P,  and  a small  increase 
in  the  other  case. 

- the  dynamic  switching  characteristics  at  a given  Reynolds 
number  are  more  dependent  on  the  output  condition  (open  or 
closed)  in  the  case  of  the  shorter  splitter  distance  than  in 
the  other  case. 

- the  best  configuration  appears  to  be  that  with  D/vv^  = 16,4, 
operated  at  moderately  low  Reynolds  numbers  (around  900,  a 
decrease  down  to  650  being  acceptable  but  not  an  increase  up 

to  1200),  and  with  control  channels  open  (or  vented  ).  /*  value  of 
about  1/400  for  the  minimum  ratio  of  control  to  supply  momentum 
flux  is  obtained.  For  a value  of  1/20  of  this  ratio,  a switching 
time  of  2.5  seconds  is  obtained.  For  a practical  device,  about 
10  times  smaller  than  the  model  and  operated  with  air  rather  than 
water,  this  switching  time,  proportional  to  w^/v,  scales  down 
to  1.7  ms.  This  date  also  indicates  that,  assuminn  a reasonable 
pressure  recovery  of  40%,  a maximum  momentum  gain  of  160  could 
be  reached. 

3 . 5 Design  and  scaling  of  the  laminar  bistable  amplif i e r 

A practical  device  will  have  much  smaller  dimensions 
than  the  model  used  here  and  will  be  operated  with  air.  As  in  the 
turbulent  flow  case,  geometric  similarity  will  be  used  to  scale 
down  the  device.  Dynamic  similarity  will  be  imposed  by  operating 
the  real  device  at  the  same  Reynolds  number  as  the  model.  Thus 
decreasing  the  bistable  dimensions  will  lead  to  increasing  supply 
velocities  and  pressures,  as  well  as  power  consumption  as 
indicated  by  the  following  relations  : 
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supply  velocity 
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volume  flow 
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pressure 
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>■  ■ Is  “p’ 

Two  limits  will  be  met  in  this  miniaturization  process  : 
technological  difficulties  in  manufacturing  small  nozzles,  of  about 
0.1  mm  width,  and  the  occurrence  of  compressible  flow  phenomena 
in  the  supply  nozzle.  This  indeed  changes  the  nature  of  the  flow 
field  in  the  interaction  region  of  the  bistable.  It  leads  to  the 
limit 


w^  > Re 


V 

rrfT 


where  C is  the  sonic  velocity,  and  '1  the  maximum  acceptable 
Mach  number.  Assuming  M=0.6,  and  for  air  at  normal  conditions 
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w^>7.3  10  Re 


For  a Reynolds  number  of  900,  this  leads  to  the  limit 
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w^  > 0,066  mm 

This  indicates  that  at  a Reynolds  number  of  900,  compres- 
sible flow  would  occur  for  a nozzle  width  smaller  than  the  smallest 
width  (0.1  mm)  that  it  is  possible  to  manufacture.  For  this  latter 
width,  compressible  flow  would  occur  when  the  Reynolds  number  • 

exceeds  1370.  j 

t 

This  also  indicates  that  selecting  a laminar  flow  bistable  ; 


instead  of  a turbulent  one  is  a necessary  step  towards  miniatur- 
ization. 


t 
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The  performance  of  a scaled  down  bistable  amplifier, 
operated  at  the  same  Reynolds  number  as  the  original  model,  is 
described  by  non-dimensional  parameters.  One  of  them  is  the 
ratio  of  a control  quantity  (flow  rate,  momentum  flux,  or 
pressure)  to  the  corresponding  supply  quantity.  Another  one  is 
the  non-dimensional  switching  tine,  which  can  be  defined  as 
Tv/w^ . The  use  of  this  parameter  instead  of  the  more  usual 
tUq/w^  is  preferred  here,  because  it  directly  relates,  for  a 
given  fluid,  switching  tire  to  geometrical  scale. 

As  an  example  the  following  table  gives  the  charac- 
teristics of  a scaled  down  laminar  bistable  amplifier.  These 
have  been  evaluated  from  the  results  of  the  tests  conducted 
with  water  on  the  bistable  model  yielding  the  best  performances, 
having  w^  = 2.5  mm  and  operated  at  a Reynolds  number  of  900. 


Nozzle  width 
(mm) 

Air  flow 
consumption 
(It/h) 

Supply 

pressure 

(kM/m2) 

Supply 
powe  r 
( mW ) 

Switching  time 
for  J /J  =1/20 
(fils)' 

0.5 

48.5 

0.45 

6 

6.7 

0.3 

29  . 1 

1.24 

10 

2.4 

0.2 

19.4 

2.79 

15 

1 . 1 

0.1 

9.7 

11.16 

30 

. 27 

A 

silhouette  of 

the  1 ami nar 

bistable 

ampl i f ier  here 

proposed  is  given  in  Fig. 48. 


3 . 6 Conclusions  and  recommendations 

A wall  attachment  bistable  amplifier  operating  in 
the  laminar  flow  regime  was  designed. 

Its  performances  compare  favourably  with  those  of  a 
turbulent  flow  bistable  amplifier.  They  are  in  the  same  range  as 
those  of  the  laminar  bistable  designed  by  Drzeviecki  at  HDL  (Ref. 12). 
The  VKI  design,  however,  is  operated  at  a larger  supply  Reynolds 


/ 


I 
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number  than  the  HDL  design.  It  has  then  a larger  power  consumption, 
for  the  same  cross  section  of  supply  nozzle  (0.5  mm^ ) , but  also 
a greater  steady  state  gain. 

Extreme  care  has  to  be  taken  in  the  fabrication  process 
to  correctly  reproduce  the  blunt  shape  of  the  splitter. 
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4.  CHEMICAL  ETCHING  TECHNIQUE  FOR  FABRICATING 
FLUIDIC  DEVICES 


4.1  Introduction 


A well-known  method  for  fabricating  devices  is  the 
chemical  etching  of  several  thin  metal  foils,  which  are  then 
piled  up  together  (Ref. 16).  A bistable  amplifier  constructed 
at  the  VKI  using  this  technique  is  shown  in  Fig. 49. 

This  process  consists  of  covering  the  metal  foil  with 
a "photoresist"  film  which  has  the  property  of  polymerizing 
under  the  action  of  UV  light.  Then  the  foil  is  exposed  to  UV 
light  through  a mask,  reproducing  the  pattern  to  be  etched. 

An  organic  solvent  is  used  to  eliminate  the  non  polymerized 
areas  of  the  photoresist  film,  and  the  metal  below  these  is 
etched  off  by  attack  of  a chemical  agent  (acid  or  alkali). 

This  process  is  reasonably  well-known  (Ref. 17).  Good 
photoresists  are  commercially  available.  Specialized  films  carry 
out  excellent  etchings  (Ref. 18),  at  the  expense  of  moderately 
large  investments. 

The  purpose  of  this  study  was  to  assess  the  possibility 
of  carrying  out  in-house  etchings  of  reasonably  good  quality, 
with  minimal  investment.  For  this  purpose,  several  simplifications 
and  modifications  to  the  processes  described  in  Ref. 17,  and  in 
particular  to  the  mask  preparation,  to  the  photoresist  treatment 
and  to  the  final  etching,  had  to  be  made. 

Etching  tests  were  then  made,  with  the  purpose  of 
determining  the  conditions  for  getting  good  quality  products  t 

(i.e.  with  minimal  and  constant  overetch)  in  the  minimal  time,  and  jj 

therefore  with  high  etch  rates.  I 
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4 . 2 Description  of  the  process  used 

4.2.1  Mask_p reparation 

As  usual,  the  mask  was  obtained  by  photographic  reduction 
of  a larger  scale  drawing  made  on  Mylar  film  or  on  any  other 
drafting  film  dimensionally  stable.  A coordi natograph  should 
normally  be  used  for  preparing  the  large  scale  drawing.  Because 
this  is  an  expensive  piece  of  equipment,  a normal  drawing  of  the 
mask  contour  was  made  with  pencil  on  a matte  drafting  film.  Black 
adhesive  tape  of  calibrated  width  was  then  used  to  manually 
fill  in  the  light  absorbing  areas  of  the  mask.  This  tape  is 
normally  available  from  the  manufacturers  of  materials  for 
electronic  printed  circuits. 

The  large  scale  artwork  was  then  photographically 
reduced  from  8:1  to  32:1,  in  several  steps.  The  VKI  photographic 
lab  was  already  equipped  for  this  kind  of  work.  The  final  mask 
was  made  of  two  identical  negative  films,  exactly  superposed 
and  fixed  one  to  the  other  along,  one  of  their  edges,  using 
adhesive  tape.  Reference  crosses  were  made  on  the  original 
drawing  to  allow  a correct  assembly  to  the  two  films  under 
microscope  observation. 

As  an  example,  a mask  used  for  the  fabrication  of 
turbulent  bistable  amplifiers  is  reproduced  in  real  dimensions 
in  Fig.  50 . 

4.2.2  Photoresist  treatment 

j 

Copper  is  the  easiest  material  to  be  chemically  etched, 
and  was  thus  selected  for  the  present  application.  As  a consequence, 
the  Kodak  liquid  photoresist  KPR  was  also  selected.  It  was 
deposited  on  a 0.1  mm  thick  copper  plate,  using  the  "dip-coating" 
technique.  The  plate  had  to  be  thoroughly  cleaned  before  applying 
this  coating.  Cleaning  was  made  with  a usual  house-hold  product 
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for  brass  and  copper.  Scratches  had  to  be  absolutely  avoided. 

A 10  seconds  dip  in  a 10  % HCl  solution  was  used  to  complete 
the  cl eani ng . 

Prebaking  the  dryed  photoresist  film,  exposing  to  UV 
light,  developping  the  photoresist  image,  and  post-baking  it 
were  carried  out  as  suggested  in  Ref. 17.  Eight  high  pressure 
mercury  vapour  lamps  of  125  watts  were  used  to  expose  the  copper 
plate,  placed  in  between  the  two  negatives  of  the  mask. 

4.2.3  Etching 

A ferric  chloride  solution  is  usually  employed  to  etch 
copper.  This  is  done  in  special  machines  where  air-less  sprays 
of  the  etching  solution  are  directed  against  the  metal  plate. 
Movement  of  the  spray  nozzles  to  uniformly  cover  the  plate  is  pro- 
vided, as  well  as  continuous  filtration  and  heating  at  constant 
temperature  of  the  etchant.  Because  the  latter  is  very  corrosive, 
all  parts  in  contact  with  it  must  be  made  of  ceramic,  plastic, 
glass,  or  titanium.  These  machines,  well  adapted  to  moderately 
large  production  rates,  are  rather  expensive. 

Another  approach  requiring  minimal  investment,  was 
adopted  at  the  VKI.  The  copper  plates  were  etched  by  simple 
immersion  in  a ferric  chloride  solution  at  ambient  temperature. 
Still  etching  baths  were  first  used.  Agitating  the  bath  by  using 
a low  speed  rotating  mixer  or  by  tilting  periodically  the  etchant 
container  was  also  tried. 

4 . 3 Results 

Straight  channels  were  etched  in  to  0.1  mm  thick  copper 
plates,  immerged  in  still  or  agitated  baths  of  ferric  chloride 
solution. 
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The  rate  of  copper  removal,  expressed  in  u/min,  and 
the  overetch,  expressed  in  u and  defined  by  comparison  of  the 
etched  plate  with  the  original  mask,  were  measured. 

These  are  given  in  Fig. 51  for  a still  bath,  as  function 
of  the  depth  below  the  liquid  free  surface.  Inspection  of  this 
figure  shows  that  : 

- an  increase  in  Fe  CI3  concentration  causes  a decrease  in  over- 
etch. Also  a decrease  in  etch  rate  is  obtained  for  concentration 
greater  than  40°B. 

- Overetch  and  etch  rate  do  not  vary  with  the  depth  of  the  copper 
plate  below  the  free  surface  level,  except  in  the  first  3 cm 
below  free  surface  level,  where  etchina  is  more  intense  and 

in  the  3 to  4 cm  above  the  bottom  of  the  Fe  CI3  container,  where 
etching  is  slow.  This  is  probably  due  to  accumulation  of  the 
by  products  of  the  chemical  reaction. 

- The  over-etch  depends  on  the  width  of  the  etched  channel. 

- A fresh  etchant  yields  higher  etch  rates  than  a used  one. 

Etching  rates  are  halved  when  approximately  2 of  copper 
per  cm^  of  Fe  CI3  solution  have  been  etched. 

- An  etch  rate  of  Ip/min  is  unacceptably  low  : it  means  that  it 
would  take  1 h 40'  to  etch  a 0.1  nm  thick  plate. 

- An  overetch  varying  from  about  PO  to  140  p is  obtained.  This 
figure,  of  the  same  order  as  the  plate  thickness,  is  unacceptably 
large. 


Still  etching  is  therefore  unsatisfactory,  and  tests 
were  thus  carried  out  by  agitating  the  etching  bath  with  a mixer. 
Two  kinds  of  mixer  were  used  : a simple  one,  L shaped,  made  from 
a bent  glass  rod  of  5 mm  diameter,  and  a second  one,  made  of  a 
rectangular  40  mm  x 100  mm  plate  glued  on  a straight  glass  rod. 
These  two  mixers  were  rotated  around  a vertical  axis  at  a speed  of 
670  RPM  and  540  RPM  respectively.  The  etching  solution  was  placed 
in  a cylindrical  container,  having  its  axis  coinciding  with  that 
of  the  mixer.  The  plates  to  be  etched  were  placed  either  parallel 
or  perpendicular  to  the  motion  of  the  etching  fluid. 
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Results  are  given  in  Fig. 52.  It  may  be  seen  that  the 
motion  of  the  etching  fluid  has  a very  strong  influence  on  the 
chemical  reaction.  Furthermore,  it  is  essential  to  have  almost 
constant  etch  rates  in  all  the  bath,  otherwise  non-uniform  and 
excessive  overetches  will  be  obtained  (curve  7).  The  trend  of 
curve  A can  be  explained  considering  that  the  L shaped  mixer 
mainly  agitates  the  lower  part  of  the  bath. 

Curve  C represents  the  best  etching  procedure.  It 
means  that  the  fluid  must  impinge  perpendicularly  onto  the 
plate.  Again  the  increase  in  overetch  towards  the  bottom  is  due 
to  the  L shaped  mixer. 

Curves  B and  D demonstrate  that  with  a mixer  which 
should  induce  near  the  plate  fluio  velocities  approximately 
constant,  reasonably  constant  etch  rates  are  also  obtained.  They 
are  however  much  lower  than  in  the  previous  cases  A and  C.  This  is 
parity  due  to  the  lower  rotational  speed  of  the  mixer,  and  partly 
to  the  higher  Fe  Cl  3 Beaume  density.  The  latter  alone  causes 
(Fig. 51)  a decrease  of  33«  in  etching  rate. 

The  following  conclusions  can  be  drawn  : 

- the  etching  bath  must  he  agitated, 

- this  agitation  must  be  uniform, 

- positioning  the  plate  parallel  or  perpendicular  to  the  motion 
of  the  fluid  has  no  effect  on  the  minimum  overetch, 

- positioning  the  plate  perpendicular  to  the  motion  of  the  fluid 
yields  higher  etch  rates. 

As  a result  of  this  study,  the  etching  apparatus  shown 
in  Fig. 53  was  constructed.  The  periodic  tilting  of  the  etching 
bath  provides  uniform  agitation,  with  velocities  parallel  to  the 
place  to  be  etched.  This  does  not  minimize  etctiing  time,  but  is 
the  siniplestand  more  economical  apparatus.  Tilting  frequency  is 
about  0.5  Hz  . 
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4.4.  Products 

The  bistable  silhouette  shown  in  Fig. 31  was  chemically 
etched  on  0.1  mm  thick  copper  plates,  by  using  the  process 
previously  described.  Several  etchings  were  made  with  nominal 
nozzle  widths  of  1.0  mm,  0.5  mm  and  0.25  mm.  A photographic 
enlargement  of  this  last  etching  obtained  using  the  mask  shown 
in  Fig. 50  i'.  reproduced  in  Fig. 54.  The  width  of  the  rectangular 
plate  is  8.90  mm  in  reality,  and  therefore  magnification  of 
Fig. 52  is  9.89.  Thus  real  width  of  supply  and  control  nozzles 
can  be  measured.  They  are  of  0.27  mm  and  0.30  mm  respectively, 
instead  of  0.25  mm  nominal.  The  difference  between  the  width  of 
the  supply  and  control  channels  is  probably  due  to  the  fact 
that  the  velocity  in  the  etching  bath  is  directed  parallel  to  the 
supply  nozzle  and  perpendicular  to  the  control  ones.  Over-etch 
is  therefore  of  10  u for  the  supply  channel  and  25  u for  the 
control  channels.  Etching  time  was  50  minutes,  and  therefore 
the  etch  rate  was  about  2 y/min.  These  figures  compare  favourably 
with  the  results  of  the  preliminary  etching  tests  described  in 
the  previous  section. 

Fig. 54  represents  a sample  out  of  a batch  of  80  units. 
14  units  in  this  batch  had  to  be  rejected  because  they  showed 
visible  defects.  The  latter  were  mainly  incomplete  etching  of 
channels.  On  the  other  hand  when  this  study  was  undertaken,  the 
defects  were  mainly  localized  over  etch,  or  pin  holes  caused 
by  bad  photoresist  adherence  to  the  plate  resulting  from  poor 
plate  cleaning  prior  to  photoresist  application.  Results 
presently  obtained  show  therefore  that  the  steps  taken  to  improve 
plate  cleaning  were  essential  and  effective. 


4.5  Conclusions  and  recommendations 


In  house  etching  of  copper  plates  with  a minimal 
amount  of  investment  proved  to  yield  products  of  reasonable 
quality.  An  improvement  to  the  process  could  be  made  by  providing 
a continuous  change  in  the  direction  of  the  velocities  in  the 
agitated  etching  bath. 
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5.  TESTING  OF  BISTABLE  AMPLIFIERS  CONSTRUCTED  BY 
CHEMICAL  ETCHING 


5 . 1 Description  of  tests 

The  purpose  of  this  study  was  to  check  the  performance 
of  the  turbulent  bistable  amplifier  previously  designed  (Fig. 31) 
and  constructed  using  a chemical  etching  technique. 

Tests  were  conducted  on  a bistable  having  a nominal 
width  of  the  supply  nozzle  of  1 mm.  The  supply,  output  and  input 
characteristics  were  determined  in  steady  state  operation. 


Static  pressures  were  measured  in  the  4 mm  internal 
diameter  plastic  tubing  connected  to  the  supply , control  and 
output  ports  of  the  bistable  amplifier.  Supply  and  control  flow 
rates  were  measured  using  two  rotameters.  Output  flow  rates  were 
obtained  from  the  measurement  of  total  and  static  pressures 
in  a portion  of  rigid  4 mm  internal  diameter  tubing  connected 
to  the  output  ports.  For  this  purpose,  a miniature  pitot  probe 
was  mounted  in  the  centerline  of  the  output  tubing. 


5.2  Test  results 


Fig.  55  shows  the  measured  supply  characteristics, 
relating  supply  flow  rate  to  supply  pressure.  The  curve  is 
dependent  on  the  configuration  of  control  channels  (open  or 
closed).  Furthermore  the  discharge  coefficient  is  rather  low, 
about  0.77  and  0.86  with  controls  respectively  open  or  closed. 

Fig. 56  shows  the  two  characteristics  measured  at  the 
two  outputs.  Again  they  depend  on  control  configuration  (open  or 
closed).  Furthermore  the  curves  for  the  two  outputs  do  not 
coincide.  This  indicates  that  the  bistable  is  not  exactly 
symmetrical,  and  therefore  that  appropriate  steps  have  to  be 
taken  to  improve  the  symmetry  of  the  etched  plates. 


The  maximum  pressure  recoveries  (0.41  with  controls 
open,  0.37  with  controls  closed)  are  smallc"  than  expected 
(0.47).  Scaling  effects  or  inaccuracies  in  the  etching  process 
account  for  these  differences. 

The  maximum  flow  recoveries  are  also  smaller  than 
expected.  This  is  due  to  a restriction  to  the  flow  cross 
sectional  area,  existing  between  the  output  duct  (P  mm^)  and 
the  outlet  tubing  (3.14  mm^ ) . 

Fig.  57  indicates  a weak  variation  of  the  maximum 
pressure  but  not  of  the  maximum  flow  recovery  with  the  supply 
pressure.  The  range  of  supply  pressures  yielding  bistable  oper- 
ation was  found  to  be  limited  (from  about  2 to  16  kN/m^). 

Fig.  58  gives  the  input  characteristics,  relating 
flow  and  pressure  at  each  control  ports.  The  characteristics 
are  dependent  on  the  output  configuration  (open  or  blocked). 
Furthermore  it  can  be  seen  that  again  the  device  is  not  exactly 
symmetrical.  Switching  of  the  bistable  occurs  for  control  flow 
rates  higher  than  expected. 

A variation  in  the  supply  pressure  between  2 to  16 
causes  almost  no  variation  on  the  ratio  of  the  control  to  supply 
flow  at  switching,  fin  the  other  hand,  the  effect  of  supply 
pressure  on  the  control  pressure  at  switching  is  illustrated 
in  Fig. 52.  Again  a strong  non  symmetricity  is  evident. 


After  these  tests  the  bistable  was  disassembled  and 
the  etched  plates  examined. 

No  apparent  asymmetry  was  noticed.  On  the  other  hand,  because  of 
the  overetch,  dimensions  were  slightly  altered  with  respect  to 
the  nominal  ones.  Indeed  the  bistable  tested  was  etched  in  a 
still  Fe  CI3  bath,  and  as  explained  in  previous  section,  this 
does  not  minimize  overetch.  Thus  nozzle  width  was  1.02  mm  instead 
of  1 mm.  This  causes  a reduction  in  the  non-dimensional  setback 
a/w^  and  length  L^^/vi  of  the  attachment  wall.  These  etching 
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inaccuracies  can  account  for  the  poor  agreement  between  the 
results  of  the  present  tests  and  of  the  previous  carried  out 
on  a large  scale  bistable  model. 

5 . 3 Conclusions  and  recommendations 

The  bistable  tested  was  not  symmetrical,  and  its 
performances  were  lower  than  expected. 

Symmetricity  could  be  obtained  by  alternating  normally 
mounted  and  inverted  etched  plates. 

On  the  other  hand  the  over  etch  has  to  be  compensated 
by  a corresponding  oversizing  of  the  mask  used  to  expose  the 
photoresist. 
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APPENDIX  A : SYSTEM  OF  EQUATIONS  DESCRIBING  THE 
JET  TO  WALL  ATTACHMENT 

A- 1 The  closed  control  case 

A- 1-1  At_an_arbitrdr^_flow_section 

The  eleven  equations  listed  in  section  1-3. 5.1  and 
derived  in  Ref. 9 are  given  here  below.  Equations  6, 7, 12b, 12c, 

14  erroneously  written  in  Ref. 9 have  here  been  suitably  corrected. 

Spread  of  the  outer  shear  layer  of  the  jet 

It  is  assumed,  as  for  a free  jet,  that 

b = E .X,  with  E = 0,274  (4) 

e e e ' ' 

Spread  of  the  inner  shear  layer  of  the  jet 

It  is  assumed,  as  for  a free  jet,  that 

b.  = E.  . X,  with  E.  = 0,274  (5a) 

if  > 0,  i.e.  if  the  jet  boundary  does  not  touch  the  wall. 
Otherwise,  Eq.(5a)  is  replaced  by  the  identity 

l.  = 0 (5b) 

and  b^  is  given  by  the  available  space  between  the  jet  and  the 
attachment  wall  (see  Eq.(7)). 

Geometrical  relation  for  the  spacing  between  the  outer  jet  edge 
and  the  wal 1 

from  Fi g . ( 23 ) 

(R  + be  + be  + tg)^  cos(a  + ^^)  = (R  + be  + bg  + cos(a  + ^^_^) 

(6) 


r 
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Here  suffixes  N-1  and  N refer  to  two  adjacent  flow  sections, 
and  4>  is  the  angle  between  the  nozzle  exit  plane  and  any  flow 
section. 

Geometrical  relation  for  the  spacing  between  inner  jet  edge 
and  attachment  wal  1 


Again  from  Fig. 23 

(R  - b^  - b.  - cos(a  - = (R  - b^  - b.  - cos(a  - 


(7) 

Eqs  6 and  7 are  here  written  in  a simplified  form  with  respect  to 
that  reported  in  Ref. 9,  where  an  unnecessary  distinction  between 
the  case  R > 0 and  R < 0 as  well  as  two  errors  in  signs  in  the 
R < 0 case  were  made.  This  had  no  consequence  on  the  results, 
because  R was  always  positive. 


Continuity  equation  for  the  outer  half  jet 


Assuming  that  the  jet  centerline  is  a streamline,  the 

mass  conservation  equation  can  be  written  as 

b +b 
c e e 


^s,e  = ^ 


u dy  ] 


where  ^ = R ^(e 

s ,e  max  . e ' 


W^/2R 


1) 


(8a) 


is  the  flow  rate  emitted  through  the  outer  half  of  the  supply 

nozzle.  Because  the  velocity  profile  is  not  symmetric,  this 

quantity  differs  from  the  analogous  one  (Q  .)  relative  to  the 

s , 1 

inner  half  of  the  jet. 


By  using  Eq.(l)  and  Eq.(3)  for  u,  the  integrals  can 
be  calculated.  This  leads  to  the  expression  : 
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- U b 

max.e  e 


R 


[ 0,45  + ^ (e  ^ - 1)]  + Q 


s ,e 


rme 


l + 0.55  b 
e e 


(8) 


Continuity  equation  for  the  inner  half  jet 


The  mass  conservation  condition  is  written  as 
0 


u dy 


-V.'  /2R 

where  = R “max.i<*  ' " > 

is  the  flow  rate  emitted  through  the  inner  half  of  the  supply 
nozzle . 


This  leads,  with  Eq.(l)  and  Eq.(3),  to  the  expression 


-U  .b.  [0.45+^(l- 
max.i  1 b,  ' 


•b./R 


)1  + Q 

^s,i 


rmi 


l.  + 0.55  b. 


(9b) 


Bernoulli  equation  for  the  external  flow 

Because  the  flow  in  the  reg’jn  outside  the  curved 
outer  boundary  of  the  jet  is  irrotational  and  cones  from  the 
still  atmosphere,  the  relative  pressure  in  such  region  is  given 
by 

Pe  = • f ^ 

where  is  the  component,  normal  to  the  flow  section,  of 

the  velocity  which  is  directed  parallel  to  the  outer  wall. 
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Additional  relation  between 

This  expression  is  obtained  by  equation  the  value 
of  Pg  given  by  Eq.(lO)  to  the  following  one  : 

b +b 
c e 

Pe  ■ Pt  ■ I "Lx.e  * I 1 

b 

c 

obtained  by  using  : 

- Bernoulli's  equation  to  predict  the  static  pressure  at  the 
potential  core  edge,  knowing  the  total  pressure  p^  of  the 
supply  f 1 ow . 

- a momentum  balance  to  predict  the  static  pressure  variation 
across  the  outer  shear  layer  of  the  jet. 

After  evaluating  the  integral  using  Eq.(l)  and  Eq.(3), 

one  gets 


e U 

2 ma  X . e 


L'2 
r me 


Pt  " 2 


COS^  ( a + i>) 


243 


(U  - U ) + 

» mav  a 


R 1^2  7 0 ' max.e  rme 


+ 0,9(U  - U ) b'  + U 

' max.e  rrne'  rme  rme 


(11) 


Momentum  balance  across  the  jet 


Consider,  in  the  inner  half  of  the  curved  jet,  the 
streamline  dividing  the  supply  flow  from  the  flow  that  has  been 
entrained  in  the  jet.  It  originates  obviously  at  the  inner  nozzle 
edge  and  terminates  at  the  attachment  point,  and  therefore  bounds, 
in  the  closed  control  case,  the  recirculating  bubble.  The  pressure 
p^  at  the  bubble  boundary  can  be  obtained  from  a momentum  balance 
across  the  jet 


Pi 


£ 

R 


dy 

-(b,^b^) 


(12a) 


where  the  Y-coordi nate  of  the  dividing  streamline,  defined  as 

- (b  + bjl  can  be  determined  from  a mass  conservation  condition 
' c d ' 


-(b.-^bd) 


(12b 


Again,  using  Eq . ( 1 ) and  Eq.(3)  to  integrate  these  expressions, 
one  gets,  from  the  momentum  balance  : 


Pe  • Pi  =Pir 


243 


™(U  - U )2  + 0,9(U 

770^  max.e  rme'  ' 


- U )U  + U 
max.e  rme'  rme 


r,  2 2 

+ ^ (U  . - U ) 

2 ' max .1  max.e 


+ P 


b. 

1 


2.5  8 55 

1.6n^  + 1.5n^  ‘ IT  ^d 


2.5 


+ n , ) + 2{U  • - U -)b  -(ij  ■ 0.4nj  + 0.25n.) 

7 d'  ' max.1  rmv  rmi'  'd  d d' 


+ 


U2  . 
rmi 


(12) 


From  the  continuity  equation  defining  n^  = b^/b^.  (Eq.(12b)) 
one  gets 


R U 


max 


.i('  - 


-br/R 
e ) 


+ U 


rmi 


b . 
1 


(U 


ma  X . 1 


U . )b.  . 
rmi  ' 1 


. (n^  - 0.8n^  + 0.25n^) 


(12< 


where  Q . is  expressed  by  Eq.(9b). 
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Momentum  conservation  in  the  jet  along  X-direction 

This  condition,  derived  in  Ref. 9 by  integrating  the 
Navi er-Stokes  equation  first  along  Y from  the  inner  to  the 
outer  jet  boundary  - and  then  along  X,  can  also  be  directly 
written  in  integral  form  as 


b +b 
e c 


W /2 
s 


(p  + pu^ ) dy  - 


(p  + pu^)  dy 


•b . -b 
1 c 


-W^/2 


x = o i 


d(b  + b ) 
' e c 


(p  + pU  )+pU 
' ^e  rme ' rme  ' e 


dx 


d(b.  + b^) 


dx 


(p.  +pU  .)  +pU  . V. 
' rmi  ' rmi  i 


dx  = 0 


where  v and  v^  , transversal  velocity  components  at  the  outer  and 
inner  jet  boundaries  are  given  by 


V = U tgf't’  + a) 
e rme  ' ' 


V.  =U  ^tg($-a) 
1 rmi  ' 


and  where  the  function  p(y),  static  pressure  in  the  jet,  is 
given  by  a momentum  balance  in  the  Y-direction 


b +b 
e c 


Pe  ■ = R J 


u^  dy , 


p^  being  given  by  Eq.(lO) 
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It  is  shown  in  Ref. 9 that  the  third  integral  in  the 

X-momentum  equation  as  well  as  the  term  pU  .v.  in  the  fourth 

rmi  1 

integral,  are  one  to  two  orders  of  magnitude  smaller  than  the 
second  integral  in  the  same  equation,  and  can  thus  be  neglected. 

Using  Eq . ( 1 ) and  Eq.(2)  to  evaluate  the  integrals, 
one  can  thus  get  the  following  relation 


pb. 


+ pb^. 


243 


(U 


770  ' max.e  rme 


243 


- U ) + 0.9  (U  - U ) U + U 

max.e  rme'  rme  rme 


770  ^^max.i  " ^rmi^^  ^^max.i  ' ^rmi  ^ "^rmi  ^rmi 


R 2 2 

+ 2b^  p.  + p^  (U  . - U ) 

c'^t  4'‘max.i  max.e' 


b 2 


2 

b.  r 

" b.Pi 


18771 


(U 


2 . 18 


1 2 
* I i 


280280  ^‘'max.e  ^rme^  70  ^^max.e  * ^rme^^rme'*’  2 ^rme 
18771  ,,,  ,,  . 18 


280280 


(U  . - U .)2  + (u  . - U lU  + i U2 
' max.1  rmi'  70  max.i  rmi  ' rmi  2 rmi 


WjPt  ■ P 


4 (^Lx.i  - ^Lx.e) 


X = 0 I 


d ( b . + b ) 

‘'’1  * “rm(>  ■ “ 


(13) 


where  the  last  integral  can  be  approximated  as  : 

^ d ( b . + b ) Nr 

<"(  * --d>  ' -i-  ■ l_(P(.o  * ♦ »c.j 
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Momentum  conservation  in  the  separation  bubble 

A control  volume  is  defined,  as  bounded  by  two 
adjacent  flow  sections,  the  attachment  wall  and  the  dividing 
streamline.  The  momentum  balance  is  written  in  the  direction  of 
the  wall,  neglecting  the  wall  shear  stresses,  as  follows 


h 


[Jf  * .V  * (‘i  * ‘■i  ■ “’d’l  „.j  - “) 


d , N 


p.  sin(a  - $)  dx^  = [J^  ^ ^i  ' ^d^l  N 


"d.N-1 

K 


d.N 


T COS(  - a)  dx 


^d.N-1 


(14a) 


for  the  N-th  control  volume,  and 


d.  1 


p^a^  cosa  + 


p^  sin(a  - t)  d 


(14b) 


^ Cjf.l  ^ Jr,i  ^ ^ b-i  - I’dl)]  (%  - «)  - 

'd.l 

T.  cosf't'  - a)  d X. 
j d ' ' d 

0 

for  the  first  control  volume. 

These  relation  differ  from  those  originally  established 
in  Ref. 3,  where  the  term  (t^+b^-b^)  was  uncorrectly  replaced 
by  t . . 
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In  these  relations  is  the  length  measured  along  the 
dividing  streamline,  and  the  momenta  of  the  forward  and  reverse 
flows  and  are  given  by  the  relation 


J ^ + J 
f r 


-b  -b  . 
c d 


-b  -b. 
c 1 


pu2  dy  + 


■b  -b  . 
c 1 


•b  -b  . - 1 . 
c ^ ^ 


pU2  . dy 

rmi 


:Pb. 
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2 5 


(U  . - U . - n . + 1 .6n  , 

' max  rmi'w70  d d 


4 8 55  1 7 

^ TT"d  ■ 7"d^  ^ 


+ U .(1  - n.) 
rmi ' d ' 


+ p 2 . U 

1 mil 


obtained  by  using  Eq.(3)  for  the  velocity  in  the  bubble. 

The  turbulent  shear  stress  along  the  dividing  stream- 
line, Tj  is,  according  to  Prandtl's  constant  diffusivity  theory 
(Ref. 7)  and  following  Sawyer's  assumption  for  curved  jets  (Ref. 10) 
given  by 


, 3u  1/  u V 
t = PC(—  - K2^) 


where  the  apparent  viscosity  e is  given  by 


G = T-  b.  (U  . - U .) 
E^.  1 ' max . 1 rmi  • 
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and  where  Eq.(3)  is  again  used  for  the  velocity  u. 

Furthermore,  the  two  integrals  in  Eq.(14)  are 
approximated  by  expressing  the  differential  term  as 

dX.  = (R  - b - b ,)  d4> 
d ' c a' 

and  by  assuming  p^  , , b^  and  b^  to  be  given  by  the 

arithmetic  average  of  their  values  at  section  N-1  and  N.  In 
t Ref. 9,  an  incorrect  expression,  with  the  term  b^  missing  was 

I adoptedinplaceoftheprevious  relation. 


Auxiliary  relation  for  the  potential  core  width 


This  relation  will  be  used  in  the  numerical  computation 
to  provide  a good  guess  for  the  value  of  b , required  to  initiate 
an  iteration  process  leading  to  the  exact  value  of  b^.  For  this 
first  guess,  the  potential  core  width  will  be  assumed  to  be  the 
same  as  that  of  a free  jet,  i.e. 


5,77  wj 


(13a) 


A-l-2At  the  first  flow  section 


At  the  first  flow  section  (X=0),  the  previously 
established  system  of  twelve  equations  degenerates  into  7 
identities  (deriving  respectively  from  Eqs  (4)(5)(6)(7)(13)(8) 
and  (9)  : 


= b.  = 0 


= a 


= a 


= W^/2 


U = U . = 0 
rme  rnn 

and  4 equations  only,  one  of  the  equations  of  the  original  system 
having  disappeared.  Indeed 
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Eq.(2)  remains  unchanged, 
Eq.(lO),  which  would  lead  to  p 


= 0,  is  replaced  by 


- ^ 

2 e 


where  v^,  the  entrainment  velocity  outside  the  jet,  is  given 


by  the  continuity  equation 


b +b 
c e 


a u . 

17 


which,  using  Eq,(l),  (3)  and  (4)  gives  for  X = 0 


W /2R 

''e  " ^max.e  ^ 'UTS'K 


0.45E) 


Eq.(ll)  reduces  to  Bernoulli  equation  at  outer  jet  edge  : 


*^e  " *^t  2 *^max.e 


- Eq . ( 1 2 ) reduce s to 


p -p.=^(U  .-U  ) 

*^e  2 ' max.1  max.e^ 


which  combined  with  Eq  . ( 1 1 ' ) , yields 


^max  . i “ *^t  *^i 


- Eq.(14)  loses  its  meaning  and  disappears. 

Furthermore,  it  is  convenient  to  replace  Eqs  (10',  10")  by  the 
following  expression 


obtaintd  by  combining  together  equations  (2),  (10*  & 10"),  (11'), 

(12-). 

Because  the  resulting  system  consists  of  four  equations 
in  five  unknowns,  an  arbitrary  value  has  to  be  given  to  one 
unknown,  p.  was  selected  as  the  quantity  to  be  fixed,  because 
experimental  results  show  that  its  ratio  to  p^,  total  jet  pressure, 
is  approximately  constant.  Indeed,  for  the  various  combinations 
of  attachment  wall  setback  and  angle  used  in  bistable  amplifiers, 
one  has  (Refs  4 and  9) 


The  following  values  have  been  used  for  the  calculation  : 


a a/W^  Pi^Pt 

10°  0.1  -0.165 

1 -0,150 

2 -0,115 

15°  1 -0,165 


15°  1 -0,165 
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A- 2 The  open  control  case 
A-2-1  Ini ti al _cgntrol _yol ume 

The  five  equations  listed  in  section  1-3. 5. 2.1  and 
derived  in  Ref. 9 are  given  here  below.  Eqs  (16)  and  (19) 
erroneously  written  in  Ref. 9,  have  here  been  suitably  corrected 

Momentum  conservation  in  the  separation  bubble 

This  equation,  which  replaces  Eq.(14b)  written  for 
the  closed  controls  case,  is  : 


p^a^  cosa  + 


p^  si n (a  - «)  dX^ 


‘^f.l  ^ 


-b  -b^ 
c d 


p(X,Y)dy 


•b  -b.-t. 
c 1 ^ 


r 

a - 1 ) - j cos(a  - 'J’)dX^  + ^ sina  (15) 


where  all  the  terms  are  calculated  as  in  section  A-1.1,  except 

- the  control  flow  momentum  J . which  is  derived  from  the 

C > 1 

experimental  correlation  = f(Q^)  established  in  Fig.  8 
of  Ref . 4 ; 

- the  cross  wise  integral  of  pressure  which  is  calculated  using 
the  value  of  pressure  given  by  a momentum  balance  across  the 
flow.  Eq.(21a),  with  n,  = 1 and  b = b.  is  thus  obtained. 

a 3 1 


Continuity  equation  for  the  outer  half  jet 

The  mass  conservation  condition  replacing  Eq.(8) 
written  for  the  closed  controls  case  is 


b +b  +i 
c e e 


L 


Q + Q 

s,e  ^ce 


u dy 


0 
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where  has  already  been  de  f i ned  ( Eq  . ( 8a  )).  Using  Eq.(l)  and 

Eq.(3)  to  evaluate  the  integrals,  one  gets 


Q - U l 
ce  rme  e 


0.45 


b U 

e ma  X . e 


b /R 
c 


+ 0.55 


b U 
e rme 


RU  (e 

max . e ' 


1) 


(16) 


The  second  hand  term  also  represents  the  volume  flow  entrained 
across  the  outer  jet  boundary  g-  This  quantity  is  used  in 
the  program  for  computation  of  intermediate  steps. 


Continuity  equation  for  the  inner  half  jet 


Again,  the  mass  conservation  condition  replacing  Eq.(9) 
written  for  the  closed  controls  case,  is  : 


‘SI 


+ Q 


Cl 


0 

f u dy 


1 


where  Q . has  already  been  defined  (Eq.(9a)).  Using  Eq.(l)  and 
(3)  to  evaluate  the  integrals,  one  gets 


b /R 
c 


Q . - U .8,.  = 0.45  b.U  . + 0.55  b.U  . + RU  .(e  - 1)  - Q . 

^ci  rmi  1 1 max.1  i rmi  max.i'  ' ^si 


(17) 


The  second  hand  term  also  represents  the  volume  flow  entrained 
across  the  inner  jet  boundary  Qg^^  This  quantity  is  used  in 
the  program  for  t h.>  computation  of  intermediate  steps. 
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Pressure  loss  in  the  outer  control  channel 


By  neglecting  the  dynamic  head,  the  relative  pressure 
outside  of  the  curved  outer  boundary  of  the  jet  is  given  by 
the  expression  : 


Pe  = - Ap(Q^g) 


(18a) 


where  Ap  = ^(Q^-g)  is  the  pressure  loss  in  the  outer  control 
channel,  associated  with  the  entrainment  of  the  flow  rate 


In  Ref. 9 the  following  relation,  derived  from  a 
numerical  fit  of  experimental  data  is  used  : 


Ap  = 0.223  V2 
^ ce 


0.563.  10*2  V‘*  + 0 . 14810"  ^ V® 

ce  ce 


0.206.10"^  V8 
ce 


where  = Q^g/Wg  (18b) 

Pressure  loss  in  the  inner  control  channel 


The  relative  pressure  outside  of  the  curved  inner 
boundary  of  the  jet  is  assumed  in  Ref. 9 to  be  given  by  the 
express i on 


2 

Urmi 

Pi  + I = - Ap(Q^i)  - AP^^^  (19) 

C0$2 ( o - $ ) 

where  the  pressure  loss  Ap  = (Qgi)  i inner  control  channel, 
is  given  as  previously,  by  Eq.(18b)  and  where  Ap^^^  is  the 
pressure  drop  associated  with  the  change  of  ^ + a in  the  flow 
direction.  In  Ref. 9,  this  term  is  approximately  evaluated  as 


A p . = 1.2 
^rot 


n/2  + 
~TIT 


Q ■ 2 

£ f Cl  \ 

2 


a 
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I 

A-2-2  At  any  flow  section  downstream  of  the  initial 
control  volume 


Out  of  the  five  equations  listed  in  the  previous 

paragraph,  two  (Eqs  18  and  19)  disappear,  as  well  as  the  two 

unknowns  Q and  Q .,  which  have  been  calculated  in  the  initial 
ce  Cl 

control  volume.  Furthermore,  Eqs  16  and  17  still  apply.  Eq.(15) 
has  to  be  modified  to  yield  Eq.(21),  as  shown  here  below.  One 
additional  equation,  replacing  Eqs  12  or  12'  for  the  previous 
cases,  appears  as  Eq.(20). 

Momentum  balance  across  the  jet 


The  pressure  p^  at  the  bubble  boundary  can  be  obtained 
from  a momentum  balance  across  the  jet. 


b +b 
c e 


_ P 


u^  dy 


-(b^.b^) 


(20) 


where  the  y-coordinate  of  the  attachment  streamline,  defined 
as  - (b^  + b ),  can  be  determined  from  a mass  conservation 

C a 

condition  similar  to  Eq.(12b)  : 


'SI 


+ Q 


Cl 


u dy  ] 


-(b  +b  ) 
' c a 


(20b) 


Using  Eq.(l)  and  Eq.(3)  to  integrate  these  expression,  one  gets, 
from  Eq.(20)  an  expression  similar  to  Eq.(12)  where  the  quantity 
n.  is  replaced  by  n,  = b /b . . Similarly,  Eq.(20b)  leads  to  the 
expression  defining  ; 


I 

j 


I 


L 


'I 
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RU  .(1  - e ) + U .b.n  + (U  .-U  .)b.(n  - 0.8n,+  0.25n,)  = 

max.i'  ' rmi  i a ' max.i  rirn  ' a a a' 


= Q . + Q . 

^S1  Cl 


(20c) 


Momentum  conservation  in  the  separation  bubble 

A control  volume  is  defined,  as  bounded  by  two  adjacent 
flow  sections,  the  attachment  wall  and  the  dividing  streamline. 
The  momentum  balance,  written  in  the  direction  of  the  wall  is, 
by  neglecting  the  wall  shear  stresses 


J , + J + 
f r 


•b  -b 
c d 


- b - b . - a . 
c ^ ^ 


p(y)  dy  cos  («>._,  - a)  + 

N-1  ^ 


p(Yd)  sin(a  - t)  dX^  = 


J . + J + 
f r 


•b  -b 
c c 


■b  -b  . -i  . 
c 1 ^ 


p(y)  dy  COS(«^  - a)  - I cos(4>|^  - a)dx 

^ X 

J ^d,N-l 


where  all  the  terms  except  the  pressure  integrals  are  calculated 
as  in  section  A-1-1. 


For  evaluating  the  integrals  in  this  expression,  the 
static  pressure  p(y)  in  a cross  section  of  the  jet  between  the 
dividing  streamline  and  the  attachment  streamline,  is  expressed 
using  the  transverse  momentum  balance 


p(y)  ■ Pi  = ^ 


-b  -b 
c a 


being  the  relative  static  pressure  on  the  attachment  streamline 
It  is  furthermore  assumed  that  the  pressure  p(y)  is,  between  the 
attachment  streamline  and  the  wall,  independent  of  y and  equal 
to  p.  . 

Then,  using  Eq . ( 1 ) and  (3)  to  evaluate  the  velocity 
integral,  the  pressure  integral  in  Eq.(21)  becomes  : 


-b  -b  . 
c d 


P(y)  dy  = p.(b.  + t.  - b^)  + (b^  - b^)^  p. 


-b  -b  . -t . 
c 1 1 


b? 

(‘^max.i 


II  % 2 1 6 3 5 , 3 5 

^mi)  i — ■ 35"  To" 


-(ti.  - nj(n 


8 2 5*  3 8 sSx 


d a'^'d  5 'd  2 'd  11  d 11 


+ 2(U  . - U . )U 

' max . 1 rmi ' i 


+ i_,^3.5  + ii 


5 rd 


rmi  2 35 


2.5  n 


■ )('3j  “ T 7 — ) -(n.  - n )n 

'd  a''d  5 d 4''  ►•mi'n  ' = > ' 


rmi  ' d a ' a 


(21a) 


where  = b^/b^ , coordinate  of  the  dividing  streamline,  is 
defined  as  in  the  closed  controls  case  (Eq.(12c))  and  where 
"a  ” ^a^^i  ’ ^oo'"dinate  of  the  attachment  streamline  is  given  by 
Eq.(20c). 
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APPENDIX  B ; THE  COMPUTER  PROGRAMS 

The  computer  programs  written  by  Pekaert  (Ref. 9) 
and  their  listings  are  given  here  below,  for  the  tv'O  cases  of 
closed  and  open  controls.  The  non  evident  meaning  of  the  variables 
used  in  the  program  is  given  in  the  list  of  symbols. 


Input  data  to  the  programs  is  : 
Geometric  parameters 
AI;  AE;  WS;  WC;  ALPHA 


)a  ramete  r s 


PTOT;  PINT  (closed  controls  case  only) 


constants 

Dll  (K,  = 9.77);  D22  (K.=7.67);  CONST  (K.,=  10) 
I 1 L. 

XN=1.5  exponent  of  velocity  profile 


initial  guess  val ucs 

QCI,  QCE  (open  controls  case  only) 

PHIM  : angular  coordinate  of  second  flow  section 
XHID  : (open  controls  case  only) 

ACC  : accuracy  criterion  for  stoppino  iterations. 


The  numerical  values  of  this  input  data  that  have 
been  used  for  the  present  calculations  are  also  given  after  the 
two  program  listings. 
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iPHIin=0;  PHIt/2  = 0;  tlR=U  XLgN1=AI;  XL1M1=AE;  flCK1=Ws72] 

UMAXI-EQ.I2^  R'-EQ.IO"',  UMAXE-EQ.2  ; QS2-EQ.8a 
QSSE-  EQ.9a  ; PPXT-^ECjiO'+lcfJDl ST  f =0 


|PHI!;1=PHIH2;  11  = 0,  K=1;  E = 0;  m=q| 


^o, 


Bcri+o.o8ws)>o 


YES  HO 

I «cTbcii  > — 


XHID  E.U2c 


IPHIN2=PHIN2  - O.SPHInI 


|Pll1^EQ.14a| 


PNi-EQ.14b 


MR  = O 


R*EQ.  12,  PHI'I  = A«W;:/  R ; WDIS".'1‘EQ.22 
AJN2  = AJ2;  Ajm=AJl;  PII  = PH1;  A19N  = A19 
BCH1=BCH;  BN1E=BWE;  BH;I=BNI; 
WPIST1=WDIST  ; XL2tn=XL3;  XLIIJ1=XLI 


NO  YE3 




GO  TO  (BO 


@ ( STOP  ) 


FLOW  CHART  OF  THE  COMPUTER  PROGRAM  FOR  THE  CLOSED  CONTROL  CASE 


FLOW  CHART  OF  BLOCK  B OF  THE  COMPUTER  PROGRAM 
FOR  THE  CLOSED  CONTROL  CASE 


XL2N1  « 0 


QESTE-EQ.I^S 


UMAXN-'J?UX£ 

UMAXN 


DAF? 


UMAXI-EQ.  2 
Q£I;t:*Eh.  17 


FL  'W  CHAHT  OF  bLO.'K  H FOP 
FOR  01  Ffi  COKTROl 


HE  Ci MPUTEB  FROOPAM 
HANNEL  ITERATION 
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STREAMLINE  PATTERN  IN  A BISTABLE  AMPLIFIER  WITH  CLOSED  CONTROL  PORTS 


BISTABLE  AMPLIFIER  WITH  OPEN  CONTROL  PORTS 
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FIG  6 b VELOCITY  PROFILES  FOR  AR=8,  A/Ws  Lw/Ws  =17 


CONTROL  CHANNELS  CLOSED 


FIG  7a  VELOCITY  PROFILES  FOR  AR  r8 , « = 10°  A/W^  = 1,  L^/Wg  = 17 
CONTROL  CHANNELS  OPEN 


FIG  8b  VELOCITY  PROFILES  FOR  AR  = 8, cY=10°  A/Wg  = 2,  Lw /Wj  = 17 
CONTROL  CHANNELS  OPEN 
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FIG  15  THE  GROWTH  OF  THE  JET  BOUNDARY  LAYER  THICKNESS  ALONG  THE  ATTACHED 
SIDE  FOR  AR^8,  Lw/Ws=17  AND  CONTROL  CHANNELS  CLOSED 
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FIG  16  THE  GROWTH  OF  THE  JET  BOUNDARY  LAYER  THICKNESS  ALONG  THE 
ATTACHED  SIDE  FOR  AR=8,  Lw/Ws=17AND  CONTROL  CHANNELS  OPEN 


Exp.  Th.  A/W 


THE  GROWTH  OF  THE  JET  BOUNDARY  LAYER  HALF  THICKNESS  ALONG  THE 
ATTACHED  SIDE  FOR  AR=8,  Lw/Ws  = 17  AND  CONTROL  CHANNELS  CLOSED 


THE  GROWTH  OF  THE  JET  BOUNDARY  LAYER  HALF  THICKNESS  ALONG  THE 
ATTACHED  SIDE  FOR  ARz8,  Lw/Ws=17  AND  CONTROL  CHANNELS  OPEN 
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THE  WALL  PRESSURE  DISTRIBUTIONS  FOR  AR  r 8,  L w /W^  = 1 7 AND  CLOSED  CONTROL  CHANN 
OPERATION 
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THE  WALL  PRESSURE  DISTRIBUTIONS  FOR  AR  =8 , Lyy /Wj  r 1 7 AND  CLOSED  CONTROL  CHANNELS 
OPERATION 


CHANNELS  OPERATION 


HE  WALL  PRESSURE  DISTRIBUTIONS  FOR  ARr8,  L^/Ws=17AND  OPEN  CONTROL 
HANNELS  OPERATION 


FIG.  29  SWITCHING  CHARACTERISTICS  OF  STRAIGHT  AND  CURVED 
WALL  BISTABLES  WITHOUT  OUTPUT  DUCTS  AND  WITHOUT 
SPLITTER 
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FIG.  30  OUTPUT  CHARACTERISTIC  OF  THE  SELECTED  DESIGN 
FOR  A TURBULENT  BISTABLE  AMPLIFIER 


FIG.  32  FIRST  CONFIGURATION  OF  THE  LAMINAR  BISTABLE 
AMPLIFIER 


FIG.  33  UNSATISFACTORY  FLOW  PATTERN  IN  THE  ' 

LAMINAR  BISTABLE  AMPLIFIER 
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FIG.  37  UNSATISFACTORY  FLOW  PATTERN  IN  THE 
LAMINAR  BISTABLE  AMPLIFIER 


FIG.  38  FLOW  PATTERN  IN  A LAMINAR  BISTABLE 
AMPLIFIER  EXHIBITING  STABILITY  UNDER 
BLOCKED  OUTPUT  CONDITION 
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Fig.  50  Typical  mask  used  for  the  fabrication 
of  turbulent  flow  bistable  amplifiers 
(scale  I:  I) 


FgCL3  solution 


FIG.  53  SKETCH  OF  THE  APPARATUS  USED  FOR 
CHEMICAL  ETCHING 


Fig.  54  Enlargement  (9,89  X)  of  an  etching 
made  in  a 0, Imm  thick  copper  plate 


OUTPUT  CHARACTERISTICS  OF  THE  CHEMICALLY  ETCHED  BISTABLE  AMPLIFIER 
t/W.:=2,  Re=  6670 
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FIG.  58  INPUT  CHARACTERISTIC  OF  THE  CHEMICALLY  ETCHED 
BISTABLE  (Ws=lmm,  t/Ws  = 2 , Re  =6670) 
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FIG.  59  VARIATION  OF  SWITCH  POINTS  WITH  SUPPLY  PRESSURE 


